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Abstract 

The  free  space  propagation  model  is  inadequate  to  predict  the  mean  path-loss  in 
ground  wireless  communication.  Also,  many  existing  propagation  channel  models  do  not 
adequately  predict  path-loss  in  rough  terrain  because  most  of  them  were  based  on  mea¬ 
surements  in  urban  areas.  Hence,  a  channel  model  that  estimates  mean  path-loss  over 
many  different  kinds  of  terrain  conditions  is  desired. 

In  this  thesis,  two  new  propagation  channel  models,  Real  Terrain  Diffraction  Model 
(RTDMOD)  and  Universal  Terrain  Channel  Model  (UTCMOD),  were  developed.  Other  ge¬ 
ometric  theory  of  diffraction  (GTD)  methods  (Epstein-Peterson  (EP)  and  Deygout  (DG)) 
agree  within  5  %  from  3  MHz  to  3GHz  with  RTDMOD.  However,  unlike  these  other 
methods,  RTDMOD  can  account  for  an  almost  unlimited  number  of  obstacles  and  has 
the  advantage  of  predicting  the  significant  diffraction  point  locations  as  well  as  the  to¬ 
tal  path-loss  for  a  given  set  of  real  terrain  data.  The  computationally  efficient  stochastic 
channel  model,  UTCMOD,  was  developed  from  numerous  Monte  Carlo  simulations  using 
RTDMOD.  UTCMOD  agrees  within  1  dB/km  over  a  broad  carrier  frequency  range  and  a 
large  terrain  height  variance  range  with  the  popular  Hata  model.  Furthermore,  UTCMOD 
is  applicable  to  many  different  kinds  of  terrain  conditions  whereas  the  Hata  model  is  most 
relevant  for  terrain  conditions  characteristic  of  large  cities  or  medium-small  cities.  Un¬ 
like  previous  models,  UTCMOD  considers  communication  link  conditions  (i.e  the  relative 
heights  of  transmitter,  receiver,  and  terrain  peaks  along  the  path)  and  provides  the  correct 
perspective  for  analyzing  the  propagation  effects.  Also  UTCMOD  can  be  used  to  optimize 
many  tactical  communication  (TAG  COMM)  and  electronic  warfare  (EW)  scenarios  for 
Low  Probability  of  Interception  (LPI)  performance. 

Since,  Low  Probability  of  Interception  (LPI)  is  important  in  both  military  and 
commercial  applications,  methods  of  analyzing  LPI  communication  have  been  developed 
throughout  the  past  decade.  However,  the  existing  LPI  formulation  is  not  sufficient  be¬ 
cause  the  previous  LPI  quality  factor,  Qlpi,  definitions  are  overly  simplistic,  too  sensitive 
to  range,  and  ignore  atmosphere  and  terrain  effects.  Fortunately,  these  problems  can 
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be  resolved  by  using  a  new  terrain  quality  factor,  Qter,  and  a  tactical  quality  factor, 
QtaCi  which  incorporates  both  system  and  environmental  parameters.  Finally,  Qter  and 
Qtac  can  be  used  to  analyze  tactical  communication  (TAC  COMM)  and  electronic  war¬ 
fare  (EW)  scenarios  in  order  to  exploit  terrain  to  create  a  communication  advantage  for 
friendly  forces.  Although  this  work  was  caste  in  a  military  perspective,  all  of  the  models, 
quality  factors,  and  concepts  are  directly  applicable  in  currently  expanding  commercial 
applications  involving  signal  propagation  in  mountainous  rural  areas. 
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Diffraction  Analysis  and  Tactical  Applications 
of  Signal  Propagation  Over  Rough  Terrain 


1.  Introduction 

Wireless  personal  communication  is  extremely  important  in  modern  warfare  to  en¬ 
hance  the  mobility  and  effectiveness  of  military  forces.  The  performance  of  wireless  systems 
is  highly  predictable  in  free  space  where  the  received  signal  power(Pr)  is  a  simple  func¬ 
tion  of  transmitted  power  {Pt),  transmitter  antenna  gain  (Gt),  receiver  antenna  gain  (Gr), 
wavelength  (A),  and  propagation  distance  (d)  from  the  source  [Sk/ar,  1988].  Specifically 


p  PtGtGA^ 

’’  (47r)2d2 


(1.1) 


However,  practical  systems  operate  in  a  non-homogeneous,  dispersive,  anisotropic, 
lossy  medium  with  random  obstructions  along  the  path.  Whereas  buildings  may  be  the 
primary  obstacles  to  radio  waves  in  a  flat  urban  environment,  hills  and  mountains  may 
be  a  major  factor  in  other  situations.  In  any  case,  obstacles  cause  additional  propagation 
loss.  In  this  work  the  effects  of  rough  terrain  are  singled  out  for  study  and  the  intervening 
medium  is  assumed  to  be  free  space. 


1.1  Background 

In  ground  communications,  the  free  space  channel  model  is  insufficient.  There  are 
many  obstacles  affecting  the  mean  path-loss.  Rough  terrain  such  as  mountains  are  an 
example  of  obstacles  along  a  ground  communication  propagation  path.  Over  a  period  of 
many  years,  a  number  of  attempts  at  predicting  the  mean  path-loss  in  the  ground  radio 
environment  have  been  made.  These  include  the  Egli  model,  the  Longley-Rice  model,  the 
Okumura  method  and  the  Hata  model  [Peterson  et  al.,  1995].  Of  these  models,  the  Hata 
model  has  often  been  used  for  analyzing  ground  communication  paths  to  predict  the  mean 
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path-loss  as  a  function  of  distance  and  antenna  height.  The  model  was  based  on  extensive 
experimental  studies  that  emphasized  the  height  of  the  receiver  point.  Hata  accounted 
for  terrain  conditions  between  the  transmitter  (Tx)  and  the  receiver  (Rx)  with  a  correc¬ 
tion  parameter  only  suitable  for  "averaged”  urban  or  rural  conditions.  Thus  a  different 
approach  is  required  if  one  wants  to  predict  path-loss  over  many  different  kinds  of  specific 
terrain  conditions. 

In  wireless  communications,  obstacles  create  multiple  propagation  paths  between  the 
Tx  and  the  Rx  which  cause  constructive  and  destructive  interference  or  what  is  commonly 
called  multi-path  interference.  It  is  especially  a  problem  in  rough  terrain  and  in  urban  envi¬ 
ronments  [  Woerner  et  al.,  1995].  During  the  past  decade,  many  researchers  have  discovered 
communication  schemes  such  as  spread  spectrum  (SS)  modulation  which  mitigate  multi- 
path  interference.  SS  provides  resistance  to  signal  interference  from  multiple  transmission 
paths  by  exploiting  the  phase  (delay)  differences  between  the  multiple  signals  arriving  over 
paths  of  variable  lengths  [Rappaport,  1996].  In  SS  the  modulated  information  is  modulated 
(spread)  a  second  time  to  expand  the  bandwidth  of  the  transmitted  signal  larger  than  the 
original  information  bandwidth  [Feher,  1995].  Thus  single  frequency  propagation  models 
must  be  extended  to  estimate  the  dispersive  effects  of  wideband  signals  such  as  SS. 

Secure  communications  that  offer  a  low  probability  of  interception  (LPI)  is  a  current 
issue  in  military  communications.  To  achieve  secure  communications,  various  technologies 
such  as  adaptive  antennas,  SS  modulation,  interference  suppression  filters,  and  error  con¬ 
trol  coding  can  be  used  to  reduce  the  probability  of  signal  interception  by  hostile  forces. 
LPI  quality  factors  (Qlpi)  [Gutman  and  Prescott^  1989]^  can  be  used  to  assess  the  per¬ 
formance  of  LPI  techniques.  Since  the  objective  of  LPI  communications  is  to  reduce,  as 
much  as  possible,  the  range  at  which  the  signal  can  be  intercepted  by  unfriendly  receivers 
while  at  the  same  time  maintaining,  or  even  increasing,  the  communication  range  of  the  in¬ 
tended  receiver,  the  terrain  between  both  friendly  and  hostile  Tx  and  Rx  can  be  exploited 
to  improve  friendly  communications  while  simultaneously  degrading  enemy  interception 
and  jamming.  As  noted  earlier,  the  terrain  between  the  Tx  and  the  Rx  alters  path- loss 
and  the  received  power,  P^.  Communication  range  depends  on  Hence,  there  is  a  need 

^LPI  quality  factors  are  derived  in  Chapter  6 
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to  incorporated  terrain  conditions  into  LPI  link  analysis  and  Qlpi-  This  analysis  will 
result  in  a  better  assessment  and  understanding  of  the  covertness  and  effectiveness  of  LPI 
communications. 


1.2  Problem  Statement  and  Objective 

Wireless  personal  communication  is  one  of  the  fastest  growing  fields  in  the  engi¬ 
neering  world  and  is  an  important  means  of  ground  communication  in  the  military.  This 
pervasive  use  of  wireless  communications  in  the  military  necessitates  continued  propaga¬ 
tion  research  and  the  development  of  accurate  propagation  models.  Since  the  free  space 
propagation  model  is  insufficient  for  analyzing  ground  communications,  many  additional 
factors  such  as  rough  mountainous  terrain  that  can  cause  additional  propagation  path-loss 
and  multi-path  interference  must  be  considered. 

Also  since  SS  is  a  commonly  used  technique  to  reduce  signal  interference  from  multi¬ 
ple  transmission  paths  and  to  improve  communication  security,  an  analysis  of  SS  propaga¬ 
tion  over  rough  terrain  is  important.  Thus  the  goal  of  this  thesis  is  to  apply  a  reasonable 
propagation  model  for  analyzing  the  mean  path-loss  of  signals  propagating  over  rough 
mountainous  terrain.  The  model  must  be  valid  over  various  frequency  ranges  and  provide 
a  valuable  tool  and  insight  for  analyzing  LPI  characteristics  of  ground  communication  sys¬ 
tems. 


1.3  Assumption 

This  research  will  be  based  on  the  following  assumptions: 

1.  Atmospheric  attenuation  will  be  neglected  because  it  can  be  readily  accounted  for 
by  the  Liebe  model  [Liehe,  1989]  [Ghordlo,  1996]. 

2.  The  system  losses  will  be  ignored  because  they  are  system  specific  and  independent 
of  the  propagation  path. 
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3.  The  terrain  conditions  are  assumed  to  be  static.  The  time  varying  case  can  be 
handled  by  repeating  the  static  calculations  at  each  instant  in  time  or  by  using  the 
terrain  statistics  as  an  input  to  the  static  model. 

4.  The  propagation  path-loss  can  be  estimated  with  a  model  that  is  a  function  of  the 
terrain  statistics.  Furthermore,  it  is  assumed  that  various  terrain  conditions  can  be 
adequately  characterized  by  standard  statistical  metrics. 

5.  The  ground  or  surface  of  the  terrain  is  assumed  to  be  a  moderately  good  conductor 
to  simplify  the  diffraction  analysis.  To  account  for  other  type  materials,  only  the 
fundamental  knife-edge  or  wedge  diffraction  calculations,  would  have  to  be  changed. 

6.  Antennas  are  assumed  to  be  omni-directional  to  simplify  the  analysis  but  antenna 
gain  can  be  easily  incorporated  into  the  link  analysis. 

1-4  Scope 

In  this  research,  the  analysis  will  be  limited  in  the  following  ways: 

1.  The  deterministic  propagation  channel  model  or  Real  Terrain  Diffraction  Model  will 
be  based  on  the  geometric  theory  of  diffraction  (GTD)  because  this  is  an  estab¬ 
lished  method  for  computing  diffracted  fields.  In  this  work,  a  ’’deterministic  model” 
means  that  the  diffracted  fields  and  hence  the  path-loss  are  computed  for  a  specific 
predetermined  terrain  profile. 

2.  Only  very  rough  terrain  conditions  (relative  to  wavelength)  between  transmitter  and 
receiver  will  be  analyzed  so  that  signal  reflections  from  planar  surfaces  can  be  ig¬ 
nored.  Furthermore,  only  forward  signal  diffraction  will  be  considered.  Although 
backward  signal  diffraction  can  be  considered  as  a  source  of  multi-path  interference, 
the  additional  time  delay  associated  with  backward  diffraction  allows  one  to  easily 
excise  it  with  an  appropriate  receiver  design. 

3.  The  stochastic  propagation  channel  model  or  Universal  Terrain  channel  Model  will 
be  developed  using  Monte  Carlo  computer  simulation  methods  [Jeruchim  et  al,  1992] 
and  path-loss  in  decibels  (dBs)  will  be  computable  from  the  terrain  roughness  statis¬ 
tics  (i.e.  mean  and  variance).  In  this  work,  a  ’’stochastic  model”  means  that  the 
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path-loss  is  estimated  from  terrain  statistics  only  and  ignores  the  actual  terrain  pro¬ 
file  details. 


4.  Interception  by  unfriendly  receivers  is  the  only  electronic  combat  situation  that  will 
be  addressed. 

1.5  Approach 

This  research  can  be  divided  into  the  following  steps: 

1.  Review  the  literature  for  effects  of  rough  terrain  on  wireless  communication. 

2.  Use  GTD  for  estimating  fields  diffracted  from  a  knife-edge. 

3.  Derive  a  method  of  applying  GTD  to  rough  terrain. 

4.  Show  that  the  knife-edge  approximation  is  reasonable  and  speeds  up  computations. 

5.  Create  a  computer  program  to  calculate  propagation  path- loss  for  a  specified  terrain 
profile  (i.e.  a  deterministic  model). 

6.  Use  Monte  Carlo  simulations  on  the  deterministic  model  to  derive  a  stochastic  prop¬ 
agation  channel  model  that  predicts  path-loss  in  decibel  (dB)  as  a  function  of  terrain 
statistics  (mean,  variance). 

7.  Modify  the  stochastic  model  to  analyze  dispersive  channel  effects  on  SS  signals. 

8.  Develop  a  Terrain  LPI  Quality  Factor  {Qter)  to  account  for  terrain  effects  LPI 
performance. 

9.  Reformulate  LPI  analysis  to  account  for  both  system  and  environmental  factors  via 
an  overall  tactical  quality  factor  (Qtac)- 

10.  Introduce  an  approach  for  simulating  tactical  communication  (TAG  COMM)  and 
electronic  warfare  (EW)  scenarios  that  accounts  for  the  effects  of  terrain  and  possibly 
exploits  terrain  to  create  a  COMM/EW  advantage. 
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1.6  Materials,  Data  and  Equipment 

No  special  provisions  are  required  for  this  research.  Only  the  existing  Sun 
SPARC  network  at  the  Air  Force  Institute  of  Technology  (AFIT)  will  be  required.  Exist¬ 
ing  software,  such  as  MATLAB  will  be  adequate  for  this  research. 

The  stochastic  channel  model  will  be  applied  to  synthetic  terrain  data  randomly 
generated  by  a  MATLAB  program.  In  addition,  real  terrain  data  of  an  area  in  Texas  will 
be  used  for  computer  simulations  that  introduce  a  new  electronic  combat  (EC)  analysis 
program  specifically  designed  to  incorporate  terrain  effects  for  electromagnetic  wave  prop¬ 
agation. 


1.7  Thesis  Organization 

This  thesis  consists  of  seven  main  chapters  and  seven  appendixes.  Chapter  11  of  this 
thesis  presents  pertinent  theory  of  signal  propagation.  The  inadequacy  of  the  free  space 
propagation  model  for  predicting  the  mean  path-loss  in  the  ground  wireless  communica¬ 
tion  are  demonstrated  with  some  examples  (e.g  multi-path  interference).  Also,  problems 
with  existing  channel  models  are  highlighted  by  reviewing  the  Hata  model,  one  of  the 
most  widely  used  propagation  models.  It  is  shown  that  the  Hata  model  inadequate  for 
many  different  kinds  of  terrain  conditions  and  is  not  valid  over  various  frequency  ranges. 
Finally,  the  chapter  ends  with  a  discussion  of  the  desired  characteristics  of  a  good  channel 
model.  Chapter  III  covers  diffraction  theory  which  is  fundamental  for  analyzing  channels 
over  rough  terrain.  The  chapter  begins  with  an  introduction  to  McNamara’s  knife-edge 
diffraction  equation  which  uses  the  geometric  theory  of  diffraction  (GTD)  to  compute  the 
total  electric  field  at  the  receiving  point  (Efotai)-  The  chapter  concludes  that  McNamara’s 
knife-edge  equation  is  a  reasonable  model  (approximation)  for  a  variety  of  wedge-like  ob¬ 
stacles.  Chapter  IV  presents  a  deterministic  channel  model  called  the  Real  Terrain  DiflFrac- 
tion  Model  (RTDMOD)  which  estimates  additional  diflfraction  path-loss  from  real  terrain 
data.  Based  on  terrain  simulations  and  detailed  diffraction  calculations  with  RTDMOD, 
a  stochastic  propagation  channel  model  called  the  Universal  Terrain  Channel  Model  is  de¬ 
veloped.  Chapter  V  presents  the  Universal  Terrain  Channel  Model  (UTCMOD)  and  it’s 
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various  applications.  UTCMOD  is  then  compared  with  the  Hata  model  to  see  how  a  model 
based  on  numerical  experiments  (UTCMOD)  compares  to  a  model  based  on  physical  ex¬ 
periments  (Hata  model).  Chapter  VI  presents  methods  for  assessing  the  LPI  performance 
of  tactical  communication  (TAG  COMM)  and  electronic  warfare  (EW)  scenarios.  By  us¬ 
ing  UTCMOD,  the  terrain  LPI  quality  factor,  Qter,  and  the  tactical  LPI  quality  factor, 
QtaCi  it  is  shown  than  terrain  can  be  exploited  to  gain  a  communication  advantage  for 
friendly  forces.  Finally,  Chapter  VII  extracts  conclusions  from  the  results  and  proposes 
recommendations  for  future  research. 

The  appendixes  were  written  for  the  convenience  of  the  reader.  Appendix  A  lists 
the  abbreviations  and  symbols  used  in  this  thesis.  For  completeness  and  to  complement 
McNamara’s  formulation,  Kraus’  derivation  for  the  total  electric  field  due  to  diffractions 
presented  in  Appendix  B.  However,  throughout  the  thesis,  McNamara’s  knife-edge  diffrac¬ 
tion  equation  is  used  to  compute  the  total  field  when  a  diffraction  point  is  present.  To 
conclude  that  the  McNamara  knife-edge  equation  is  a  reasonable  model  (approximation) 
for  a  variety  of  wedge- like  obstacles,  the  Balanis’  wedge  diffraction  equation  is  used  to  com¬ 
pare  computations.  For  the  reader’s  convenience,  the  Balanis’  wedge  diffraction  equation 
is  derived  in  Appendix  C.  Appendix  D  presents  a  new  method  that  estimates  lit  diffraction 
path-loss  (Liit)  over  multiple  subpath  obstacles  using  path  clearance  and  the  first  Fresnel 
zone  radius  as  model  parameters.  The  MATLAB  codes  used  in  this  thesis  can  be  found  on 
the  world  wide  web  at  www. afit.af.mil/Schools/EN/ENG/LABS/C3EMRT/jyc.htm.  Ap¬ 
pendix  E  contains  the  code  names  and  details  for  obtaining  the  code  online  on  the  Internet 
for  those  interested  in  follow-on  research.  Appendix  F  is  a  bibliography.  Finally,  this  the¬ 
sis  includes  a  useful  index  arranged  in  alphabetical  order  for  the  reader’s  convenience  in 
Appendix  G. 

Also  note  that  although  this  work  is  caste  in  a  military  perspective,  all  of  the  mod¬ 
els,  quality  factors,  and  concepts  are  directly  applicable  in  currently  expanding  commercial 
applications  involving  signal  propagation  in  mountainous  rural  areas. 
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II.  Pertinent  Theory  of  Wireless  Communication 

2. 1  Overview 

A  fundamental  requirement  for  wireless  communication  is  to  maintain  a  sufficiently 
low  probability  of  bit  error  (Pg)  which  is  a  function  of  the  received  signal  power  (Pr)- 
Since,  propagation  obstacles  that  cause  propagation  path-loss  and  multi-path  interference 
can  adversely  effect  Pg,  improvements  in  wireless  communication  require  a  reduction  or 
mitigation  of  propagation  path-loss  and  the  multi-path  interference. 

In  this  chapter,  pertinent  topics  in  signal  propagation,  spread  spectrum,  and  channel 
modeling  will  be  reviewed.  A  special  technique  called  the  SS  modulation  is  widely  used  in 
wireless  communication  because  of  it’s  inherent  capability  to  reduce  multi-path  interfer¬ 
ence.  Furthermore,  a  number  of  attempts  at  predicting  the  mean  path-loss  in  the  ground 
radio  environment  have  been  made.  Of  these  models,  the  Hata  model  is  often  used  for 
analyzing  ground  communication  paths. 

2.2  Signal  Propagation 

If  there  were  no  atmosphere  nor  obstacles  in  the  signal  propagation  path,  the  wireless 
communication  propagation  problem  would  be  much  less  complex.  However  since  wireless 
communication  signals  are  attenuated  by  the  earth’s  atmosphere  and  obstacles,  their  ef¬ 
fects  must  be  included  to  account  for  the  signal  loss  observed  on  practical  communication 
links  [O’Brien  et  al.,  1985].  In  addition  to  transmitter  signal  power  (Ft),  absorption  and 
scattering  of  signal  power  by  the  earth’s  atmosphere  and  obstacles  along  the  propagation 
path  are  primary  factors  that  determine  the  communication  distance.  We  call  this  signal 
power  reduction  “path-loss”. 


2-1 


In  addition,  constructive  and  destructive  interference  occurs  when  there  are  multi¬ 
ple  propagation  paths  between  the  Tx  and  the  Rx.  This  is  called  multi-path  interference. 
Although  multi-path  interference  does  not  occur  in  free  space,  practical  ground  communica¬ 
tion  propagation  paths  require  analysis  of  the  multi-path  problem.  Multi-path  interference 
is  especially  a  problem  in  rough  terrain  and  in  urban  environments  due  to  the  presence 
of  numerous  large  scattering  objects  such  as  mountains  and  buildings  [  Woerner  et  al.,1995]. 


2.2.1  Received  Power.  In  free  space,  the  received  signal  power(Pr)  depends 
on  parameters  such  as  transmitted  power,  antenna  gain,  wavelength,  and  distance  [Sklar, 
1988]. 


(47rd)2 

•  Pt  :  Transmitter  power 

•  G<  :  Transmitter  antenna  gain 

•  Gr  :  Receiver  antenna  gain 

•  A  :  Wavelength 

•  d  :  Distance  between  transmitter  and  receiver 

Furthermore,  frequency  is  inversely  proportional  to  wavelength.  That  is 


(2.1) 


where  Vp  is  the  velocity  of  propagation  determined  by  medium.  Thus  in  free  space,  radio 
waves  propagate  with  a  path-loss  characteristic  of  the  form  A{f)d~'^.  In  ground  communi¬ 
cations,  the  free  space  propagation  model  is  insufficient.  There  are  many  factors  affecting 
the  mean  path-loss.  Over  a  period  of  many  years,  a  number  of  attempts  at  predicting  the 
mean  path-loss  in  the  ground  radio  environment  have  been  made.  These  include  the  Egli 
model,  the  Longley-Rice  model,  the  Okumura  method,  and  the  Hata  model  [Peterson  et 
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al.,  1995]  which  are  discussed  in  section  2.4. 


2.2.2  Multi-Path  Interference.  In  ground  communications,  complex  prop¬ 
agation  paths  cause  multiple  signals  to  arrive  at  the  Rx  as  shown  in  Fig. 2.1.  Although 
the  signal  originates  from  only  one  transmitter,  the  car  receiver  is  subjected  to  numerous 
signals  that  propagate  through  different  paths.  Similarly,  rough  terrain  can  create  many 
indirect  paths  in  addition  to  the  direct  or  intended  path.  The  signals  over  these  multiple 
paths  interfere  with  each  others  [Parsons,  1992].  This  phenomena  is  called  “multi-path 
interference” .  If  the  receiver  detects  the  multi-path  signals,  then  the  received  signal  is  very 
different  from  the  transmitted  signal  and  this  causes  degraded  analog  signals  or  decision 
errors  in  digital  signals. 

One  way  to  mitigate  multi-path  interference  is  with  spread  spectrum  communication 
techniques.  The  next  section  reviews  spread  spectrum  communication. 


Figure  2.1  Multi-path  from  the  transmitter  to  the  receiver  [Parsons,  1992] 
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2.3  Spread  Spectrum 

In  all  communications  systems,  the  modulated  waveform  occupies  a  frequency  band¬ 
width  that  is  dependent  upon  the  modulation  method  used  and  the  data  being  sent.  The 
modulation  scheme  under  consideration  for  this  research  is  spread  spectrum  (SS)  commu¬ 
nications.  This  field  grew  out  of  research  efforts  during  World  War  II  to  provide  a  secure 
means  of  communication  in  hostile  environments  [O’Brien  et  al,  1985]. 

In  a  SS  system,  the  transmitted  bandwidth  of  the  signal  has  been  “spread”  over  a 
larger  bandwidth  than  the  original  signal  requires.  To  be  classified  as  a  SS  system,  the 
communication  modem  must  have  the  following  characteristics  [Rappaport,  1996]: 

1.  The  transmitted  signal  energy  must  occupy  a  bandwidth  which  is  larger  than  the 
information  bit  rate  (usually  much  larger)  and  which  is  approximately  independent 
of  the  information  bit  rate. 

2.  Demodulation  must  be  accomplished,  in  part,  by  correlation  of  the  received  signal 
with  a  replica  of  the  signal  used  in  the  transmitter  to  spread  the  information  signal. 

SS  has  several  applications  and  desirable  advantages  such  as  improvement  of  multi-path 
interference  rejection  capabilities,  code  division  multiple  access  (CDMA)  applications,  and 
enhanced  low  probability  of  interception  (LPI)  capabilities^  [Feher,  1995].  As  an  example, 
the  U.S.  Army  uses  modern  tactical  communication  systems  designed  with  SS  technology 
precisely  because  these  systems  allow  many  users  to  operate  in  the  same  band  without  ex¬ 
cessive  interference  [iSass,  1983].  In  this  research,  only  the  multi-path  interference  rejection 
capabilities  of  SS  is  of  interest.  Also  since  propagation  effects  are  frequency  dependent,  one 
must  consider  how  channel  dispersion  affects  SS  performance.  Therefore,  the  multi-path 
interference  rejection  capability  of  SS  and  the  bandwidths  of  a  variety  of  SS  signals  are 
discussed  briefly  in  the  next  section. 


2.3.1  Rejection  Capability  of  Multi-path  Interference.  Basically,  the  SS 
system  is  a  system  in  which  the  transmitted  signal  is  spread  over  a  frequency  much  wider 

’^LPI  capabilities  will  be  discussed  in  chapter  6 
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than  the  minimum  bandwidth  required  to  send  the  signal.  The  baseband  signal  is  spread 
out  with  a  spreading  code  sequence.  The  actual  signal  spreading  may  be  achieved  with 
one  of  three  basic  techniques.  These  include  : 

1.  Direct  Sequence  SS  (DSSS)  in  which  the  information  is  spread  by  multiplying  the 
information  with  a  high  bit  rate  pseudo-random  digital  code  sequence. 

2.  Frequency  Hopping  SS  (FHSS)  in  which  the  carrier  is  pseudo-randomly  shifted  into 
a  number  of  different  cells  of  bandwidth  equal  to  that  of  the  information  bandwidth. 
The  shifting  is  controlled  by  a  code  sequence. 

3.  Time  Hopping  SS  (THSS)  in  which  bursts  of  the  modulated  information  are  trans¬ 
mitted  at  pseudo-random  time  intervals  dictated  by  a  code  sequence^. 

As  an  example  of  how  SS  can  mitigate  multi-path  interference,  DSSS  will  now  be  described 
in  more  detail. 

DSSS,  a  popular  form  of  spread  spectrum  processing  is  performed  by  eflfectively 
multiplying  the  original  signal  by  a  pseudo-noise  (PN)  digital  signal.  The  PN  signal  is  a 
fast  code  that  is  many  times  the  data  rate.  They  are  called  “pseudo-noise”  because  they 
are  not  real  random  Gaussian  noise  but  are  deterministic  codes  that  appear  to  be  random. 

First  the  PN  code  is  modulated  onto  the  information  signal  using  one  of  several 
modulation  techniques  (eg.  BPSK,  QPSK,  etc  ).  Then,  a  doubly  balanced  mixer  is  used 
to  multiply  the  RF  carrier  and  PN  modulated  information  signal.  This  process  causes 
the  RF  signal  to  be  replaced  with  a  very  wide  bandwidth.  The  demodulation  process 
then  simply  consists  of  the  mixing/multiplying  of  the  same  PN  modulated  carrier  with 
the  incoming  signal.  Fig.2.2  shows  this  basic  DSSS  technique.  The  output  of  the  mixer 
in  the  receiver  is  a  signal  that  is  a  maximum  when  the  two  signals  are  exactly  equal  to 
one  another  or  are  “correlated.”  The  correlated  signal  is  then  filtered  and  sent  to  the 
demodulator.  That  is  the  correlation  property  of  SS  codes. 

In  ground  communication,  rough  terrain  causes  one  direct  path  and  any  number 
of  indirect  paths.  If  the  transmitter  and  receiver  use  the  same  spread  code,  then  the 

^THSS  is  not  addressed  in  this  research  because  THSS  systems  have  found  no  commercial  application 
to  date.  However,  the  arrival  of  cheap  random  access  memory  (RAM)  and  fast  micro-controller  chips  make 
time  hopping  a  viable  alternative  SS  technique  for  the  future. 
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Figure  2.2  Basic  SS  technique. 


receiver  is  synchronized  with  the  direct  path  signal  spread  code  and  communication  is 
possible.  However,  the  received  multi-path  signal  spread  code  will  be  out  of  phase  with 
the  transmitted  code  because  the  multi-path  received  signal  has  a  different  time  delay.  Due 
to  the  correlation  properties  of  SS  code,  the  multi-path  received  signal  will  be  weak  and 
effectively  minimizes  the  multi-path  interference  problem  [Parsons^  1992]. 

Fig.2.3  shows  the  received  power  density  of  the  two  different  cases.  In  one  case,  the 
code  is  matched;  in  the  other  case,  the  code  is  unmatched.  The  direct  path  signal  with 
the  matched  code  has  strong  received  power  in  the  receiver  (Fig. 2. 3. a).  In  contrast,  the 
multi-path  signal  with  the  unmatched  code  has  much  less  power  (Fig.2.3.b).  The  receiver 
ignores  the  multi-path  received  signal  and  detects  the  direct  path  signal.  In  this  manner, 
the  multi-path  interference  problem  can  be  mitigated  using  SS  techniques  [Rappaport, 
1991]. 


2.3.2  Spread  Spectrum  Bandwidth.  The  bandwidth  in  DSSS  systems  is 
often  taken  as  the  null-to-null  bandwidth  of  the  main  lobe  of  the  power  spectral  density 
plot.  The  null-to-null  bandwidth  of  this  lobe  is  2Rc,  where  Rc  is  the  chip  rate.  Therefore, 
the  bandwidth  of  a  DSSS  system  is  a  direct  function  of  the  chip  rate.  Fig.2.4  shows  a  basic 
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1.0, 


Code  Matched 
(Direct  Path  Signal) 


1.0, 


Additional  Time  Delay 
:Tc 


Tc  Tc  Tc  To  Tc  To  Tc  To  To  Tc 


(a)  code  matched  (b)  code  unmatched 

Figure  2.3  Power  spectra  of  the  received  signal 

DSSS  signal  and  power  spectrum  density. 

In  contrast  to  DSSS,  the  bandwidth  of  FHSS  is  determined  by  the  number  of  hop 
frequencies  because  frequency  hopping  relies  on  frequency  diversity  to  combat  interference. 
This  is  accomplished  by  a  multiple  frequency,  code  selected,  FSK  (frequency  shift  keying) 
technique.  Basically,  the  incoming  digital  stream  is  shifted  in  frequency  by  an  amount 
determined  by  a  code.  This  effectively  spreads  the  signal  power  over  a  wide  bandwidth. 
One  can  think  of  the  FHSS  transmitter  as  a  pseudo-noise  PN  code  controlled  frequency 
synthesizer.  Fig.2.5  shows  a  typical  FHSS  modem.  The  instantaneous  frequency  output  of 
the  transmitter  jumps  from  one  value  to  another  based  on  the  pseudo-random  input  from 
the  code  generator.  Varying  the  instantaneous  frequency  results  in  an  output  spectrum 
that  is  effectively  spread  over  the  range  of  frequencies  generated.  In  this  system,  the 
number  of  discrete  frequencies  determines  the  bandwidth  of  the  system.  Fig.2.6  shows  a 
pictorial  representation  of  FHSS  where  is  the  signal  bandwidth  before  spreading  (s^ 
(t)  in  Fig.2.5)  and  is  the  signal  bandwidth  after  spreading  {st  (t)  in  Fig.2.5). 

Whereas  in  a  DSSS  system,  bandwidth  is  a  direct  function  of  the  chip  rate,  in  a 
FHSS  system,  bandwidth  is  determined  by  the  number  of  discrete  frequencies  that  in  turn 
is  determined  by  SS  the  code.  Hence,  bandwidth  can  be  estimated  by  an  operator  who 
knows  the  properties  of  the  SS  code  in  use.  Also,  since  the  lowest  frequency  (/,)  and  the 
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d(t)  dataratefR) 

data  =  1/T 0it/sec) 


‘'W  dataratefRc) 

Spreading  _ 

code 


C|(t)c(l)  dataratefRc) 

=  i/rc 


Power  spectrum  density  i 
of  data -and  \ 

spreading  code-modulated  i 
carrier  i 


Figure  2.4  Basic  DSSS  signal  and  power  spectrum  density. 
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(a)  Transmitter 


Estimated 

data 


(b)  Receiver 


Figure  2.5  Basic  FHSS  modem  [Peterson  et  al.,  1995]. 


Time 


a)  Transmitted  sigani 
4-ary  FSK 


b)  Transmitted  signai 
4-ary  FSK  and  FHSS 


Figure  2.6  Pictorial  representation  of  FHSS  [Peterson  et  al.,  1995]. 


2-9 


highest  frequency  (fh)  can  be  set,  the  propagation  bandwidth  is  equal  to  fh,  -  fi-  For 
convenience,  DSSS  and  FHSS  parameters  are  summarized  below: 


1.  DSSS. 

•  bandwidth:  2i?c(Hz).  where  Rc  is  chip  rate. 

•  frequency  range:  [/<.  -  i?^,]  ~  [fc  +  Rc\-  where  /c  is  carrier  frequency. 

•  fl-  [fc  -  Rc] 

•  fh-  [fc  +  Rc] 

2.  FHSS. 

•  bandwidth:  k-  Af. where  k  is  the  number  of  discrete  frequencies  that  determined 
by  SS  code,  Af  is  hop  frequency  spacing. 

•  frequency  range:  /<,  ~  [fc  +  k-  Af  ].  where  fc  is  carrier  frequency. 

•  fl-  fc 

•  A:  [fc  +  k-  Af] 

2.3.3  Propagation  of  SS  Signals.  Due  to  the  relatively  large  bandwidth  asso¬ 
ciated  with  SS,  performance  can  be  degraded  in  dispersive  channels.  Although  SS  methods 
are  often  effective  in  reducing  multi-path  interference  problems,  the  performance  of  SS  in 
extraordinarily  complex  propagation  paths  such  as  those  characterized  by  extremely  rugged 
terrain  is  not  well  understood  [  Woener  et  al,  1995].  Rough  terrain  has  many  different  prop¬ 
agation  obstacles  adding  complexity  to  path-loss  predictions  for  wide  bandwidth  signals. 
Analyzing  SS  propagation  over  complex  terrain  is  an  important  step  towards  improving  SS 
performance  in  ground  communications.  Therefore,  previous  channel  models  that  estimate 
path-loss  are  discussed  in  the  next  section. 


2.4  Channel  Models 

The  prediction  of  path-loss  is  a  very  important  step  in  planning  an  effective  ground 
wireless  radio  system.  Hence,  much  research  focuses  on  creating  propagation  models  that 
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determine  the  relationship  between  terrestrial  paths  and  path-loss.  There  are  several  prop¬ 
agation  models  primarily  developed  for  point-to-point  mobile  communication.  These  mod¬ 
els  are  based  on  experimental  results  and  are  exemplified  by  the  Egli  model  [Egli,  1957], 
the  Longley-Rice  model  [Longley  and  Rice,  1968],  the  Okumura  method  [Okumura  et  al, 
1968],  and  the  Hata  model  [Hata,  1980],  which  are  all  statistical  propagation  models.  For 
example,  the  Egli  model  is  a  statistical  model  for  predicting  propagation  loss  in  an  urban 
or  rural  environment.  However,  it  does  not  include  diffraction  losses  caused  by  propagation 
over  irregular  terrain.  The  Longley-Rice  model  is  applicable  to  point-to-point  communi¬ 
cation  systems  in  the  frequency  range  from  40  MHz  to  100  GHZ.  However,  it  does  not 
consider  correction  factors  to  account  for  the  effects  of  buildings.  The  Okumura  method  is 
based  on  an  extensive  experimental  study  of  propagation  path-loss  predictions  in  various 
terrains  for  a  specific  frequency  range.  Okumura  implicitly  takes  into  account  the  effects  of 
building  in  urban  environments.  The  Hata  model  which  is  based  on  the  Okumura  method 
is  a  widely  used  propagation  model  that  predicts  path-loss  in  urban  areas  as  a  function  of 
the  communication  distance  and  antenna  heights  in  the  following  way  [Delisle  et  al.,  1985] 
[Jeruchim  et  al.,  1992]  [Peterson  et  al.,  1995]: 

L  =  69.55  -1-  26.16/o5f(/<,)  -  13.82log{hu)  -  a{hre,  fc)  +  [44.9  -  %.hhlog{hte)]  +  log{d)  (2.2) 

•  L  :  Path  loss  (dB) 

•  fc  :  Carrier  frequency  (MHz) 

•  hte  '.  Transmitter  Antenna  Height  (m) 

•  hre  :  Receiver  Antenna  Heigh  (m) 

•  a{hre,  fc)  •  Correction  factor  (dB).  Which  is  an  explicit  function  of  receiver  antenna 
height  and  carrier  frequency. 

•  d  :  Communication  distance  between  transmitter  and  receiver  (m) 

where  the  correction  parameter,  a,  is  used  to  account  for  variation  in  the  Rx  antenna, 
height.  Expressions  for  the  correction  parameter  for  large  cities  and  medium  to  small 
cities  where  developed  from  measured  data  and  are  given  below  for  easy  reference. 
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o^iKcfc)  =  [l-1^05(/c)  -  0.7]/ire  -  [1-56Zo5(/c)  -  0.8](dB)  for  a  medium-small  city 

0({hre,fc)  =  8.29{log{l.54:hre)f  -  for  a  large  city, /c  <  300Miif2 

(^iKe,fc)  =  3.2[/osi(11.75/ire)]^  -  4.97(dj5)  for  a  large  city, /c>  300Mii'2; 

With  the  Hata  model,  path-loss  can  be  computed  from  antenna  heights,  the  communica¬ 
tion  distance  between  Tx  and  Rx,  and  the  carrier  frequency.  Since  these  parameters  are 
usually  known,  the  Hata  model  provides  an  easy  method  for  predicting  path-loss  [Peterson 
et  al.,  1995].  However  the  Hata  model  is  not  sufficient  for  predicting  path-loss  in  rough 
terrain.  The  model  is  most  relevant  for  terrain  conditions  characteristic  of  large  cities  or 
medium-small  cities.  Hence,  the  correction  parameter  is  inadequate  for  many  different 
kinds  of  terrain  conditions  not  studied  by  Hata. 


2.5  Summary 

The  prediction  of  propagation  path-loss  is  a  central  question  in  the  planning  of  wire¬ 
less  communication  services  and  a  number  of  approaches  and  statistical  models  are  avail¬ 
able  for  the  prediction  and  calculation  of  path-loss  for  a  variety  of  conditions.  However, 
most  propagation  models  have  not  been  developed  specifically  for  application  to  a  specific 
communication  channel,  but  rather  from  a  more  general  perspective  [Delisle  et  al.,  1985]. 
In  ground  communications,  rough  terrain  presents  a  variety  of  propagation  obstacles  which 
increase  the  complexity  of  the  direct  path.  Since  communication  channels  over  rough  ter¬ 
rain  are  highly  variable,  a  channel  model  that  estimates  mean  path-loss  over  many  different 
kinds  of  terrain  conditions  is  desirable. 

In  addition,  analysis  of  SS  signal  propagation  over  rough  terrain  using  channel  mod¬ 
els  that  explicitly  account  for  terrain  statistics  is  an  important  step  towards  improving  SS 
performance  in  rough  terrain  environments. 

In  the  next  chapter,  the  effect  of  obstacles  on  the  propagation  of  elctromagnetic  fields 
will  be  discussed  and  these  results  will  form  the  basis  for  a  new  propagation  model. 
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III.  Diffraction  Theory 


3. 1  Overview 

A  radio  systems  engineer  may  be  interested  in  knowing  a  communication  system’s 
area  coverage  in  areas  that  are  highly  mountainous  or  contain  other  irregularities  that 
invalidate  the  standard  planar  earth  approximations  typically  used  for  computing  path 
attenuation.  With  the  availability  of  fast  computers,  it  is  now  possible  to  make  numer¬ 
ical  estimates  of  path-loss  for  rough  terrain  using  the  Geometric  Theory  of  Diffraction 
(GTD),  a  technique  that  evaluates  diffraction  phenomena.  In  this  thesis,  calculations  will 
be  limited  to  diffraction  from  perfect  electronical  conducting  (PEC)  bodies  to  simplify  the 
diffraction  analysis. 

Normally  signal  diffraction  occurs  when  an  obstacle  exists  in  between  the  Tx  and 
the  Rx.  If  the  obstacle  blocks  the  direct  line  of  sight  (LOS),  then  the  type  of  diffraction 
is  called  “shadow  diffraction”.  On  the  other  hand,  if  the  obstacle  is  below  LOS,  then  the 
diffraction  is  called  “lit  diffraction”.  In  this  chapter,  lit  diffraction,  a  common  cause  of 
multi-path  interference,  and  shadow  diffraction,  a  typical  contribution  to  path-loss,  will  be 
analyzed  in  different  ways. 


3.2  History  of  Diffraction  Theory 

In  1947,  Bullington  suggested  that  a  series  of  obstacles  can  be  replaced  by  a  single 
equivalent  obstacle  as  a  means  of  obtaining  the  path-loss  in  a  simple  manner  [Bullington, 
1947].  Prom  that  idea,  numerous  methods  for  computing  attenuation  due  to  diffraction 
were  developed  [Pogorzelski,  1982]. 

New  approaches  using  the  Geometric  Theory  of  Diffraction  (GTD)  have  received 
considerable  attention  in  the  past  few  years.  These  methods  provide  a  correct  diffraction 
loss  factor  and  yield  results  consistent  with  the  predictions  of  other  methods.  The  GTD, 
originated  by  Keller  [Keller,  1962]  and  extended  by  Pathak  and  Kouyoumjian  [Pathak  and 
Kouyoumjian,  1974]  has  been  combined  with  the  moment  method  for  analyzing  structures 
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which  are  small  in  terms  of  wavelength  [Burnside  et  al.,  1975]. 


3.3  McNamara  Diffraction  Equation 

3.3.1  Huygens’  Principle.  Huygens’  principle  states  that  each  point  on  a  pri¬ 
mary  wavefront  can  be  considered  to  be  a  new  source  of  a  secondary  spherical  wave  and 
that  a  secondary  wavefront  can  be  constructed  as  the  envelope  of  these  secondary  spherical 
waves,  as  suggested  in  Fig.3.1.  This  fundamental  principle  of  physical  optics  can  be  used 
to  explain  the  apparent  bending  of  radio  waves  around  obstacles,  (i.e.,  the  diffraction  of 
waves).  A  diffracted  ray  is  one  that  follows  a  path  that  cannot  be  interpreted  as  either 
reflection  or  refraction. 


Figure  3.1  Illustrating  Huygens’  principle  of  physical  optics  (point-to  wave  correspon¬ 
dence)  [Kraus,  1992]. 

3.3.2  Total  Received  Electric  Field  {Etotai)  •  Fig.3.2  illustrates  the  geometry 

for  analyzing  diffraction  from  a  perfectly  conducting  half  plane.  By  using  Huygens’  prin¬ 
ciple,  the  total  field,  Etotai,  at  an  observation  point  in  the  region  x  >  0  can  be  expressed 
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Figure  3.2  Diffraction  from  a  perfectly  conducting  half-plane  illuminated  by  a  plane  wave 
at  normal  incidence. 


in  the  following  manner  [McNamara  et  al.,  1990] 


EqC 


— jr/4)  foo 
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for  X  >  0  where 


(3.1) 

(3.2) 


X  =  distance  from  x=0  reference  plane  containing  the  obstacle  (conducting  half-plane), m 
y  =  distance  into  shadow  or  lit  region, m 
A  =  wavelength,m 
Eg  =  free  space  field  strength 
k  =  2'k/X 

7  =  2/\/^. 

^Kraus  also  derived  an  equation  for  total  field.  See  Appendix  B.  Kraus  Diffraction  Method  for  more 
details. 
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The  Fresnel  integrals  in  Eq.3.2  are  defined  as  follows  [Abramowitz  and  Stegun,  1964]: 


C{x) 

S{x) 


-KV? 

COS - du 

2 

Sin - du 

2 


Fresnel  cosine  integral 
Fresnel  sine  integral. 


If  the  obstacle  blocks  the  direct  line  of  sight  (LOS)  between  the  Tx  and  the  Rx,  only 
the  diffracted  signal  can  be  detected  because  Rx  is  in  the  shadow  region.  In  this  case,  the 
total  field  is  solely  due  to  shadow  diffraction. 

However,  if  the  obstacle  does  not  block  the  direct  LOS,  then  the  receiver  detects 
both  a  diffracted  field  (  called  Lit  diffraction)  and  the  original  incidence  field.  Thus  these 
two  cases,  shadow  and  lit  diffraction,  will  be  distinguished  throughout  this  thesis. 

For  easy  reference,  the  different  diffraction  cases  are  summarized  in  Table  3.1^. 


Diffraction 

Boundary 

Position  of 

Rx 

Direct  LOS 
between  Tx  and  Rx 

Etotai 

Shadow 

Diffraction 

Rx  in 

Shadow  region 

Blocked 
by  Obstacle 

^diffraction 

Lit 

Diffraction 

Rx  in 

Lit  region 

Not  Blocked 
by  Obstacle 

Einddence  T  Edif fraction 

Table  3.1  Diffraction  boundary 


3.3.3  X  and  Y  Values.  In  Eq.3.2,  total  field  {Etotai)  the  Rx  can  be  calculated 
with  X  (distance  from  reference  plane  containing  the  obstacle)  and  y  (distance  into  shadow 
or  lit  region).  If  the  incidence  field  (Einddence)  is  not  perpendicular  to  the  obstacle,  one 
must  redefine  the  x  and  y  values  as  shown  in  Fig.3.5  before  using  Eq.3.2  to  calculate  Etotai- 
Thus  for  an  arbitrary  incidence  angle,  Etotai  can  be  computed  as  follows: 

^The  total  field  in  the  x  <  0  region  would  be  a  superposition  of  the  incident  and  diffracted  fields  along 
with  a  component  refiected  from  the  obstacle. 
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(3.3) 


E^e-^{kx-K/4) 

V2 

where  x  and  y  are  now  given  by 

X  =  p  •  cos[7r  —  {(j)  —  (f>)\ 

=  p  •  cos\{^  —  (j)')  —  tt] 
y  =  p  •  sin[7r  —  {(f)  —  <f)')] 

=  —p  •  sin[{(f)  —  (f)')  —  tt] 


{[^-Ci-yy)]-j[--S{-jy)]} 


in  the  lit  region 

(3.4) 

in  the  shadow  region 

(3.5) 

in  the  lit  region 

(3.6) 

in  the  shadow  region. 

(3.7) 

Note  that  the  lit  and  shadow  regions  are  now  defined  in  the  following  manner 


Figure  3.3  Incidence  angle((^’)  and  Observation  angle((^). 


Lit  Region  :  7r/2  <  ((^  —  (/>')<  tt 

Shadow  Region  :  tv  <  {(f)  —  <f)')  <  2n 


(3.8) 


3.3.4  Wedge  Obstacles.  Although  Eq.3.2  gives  Etotai  observation  point 

and  is  easily  used  as  a  key  concept  for  estimating  signal  path-loss,  it  is  based  on  the 
assumption  that  the  obstacle  is  a  knife-edge.  However,  normally  the  obstacle  is  not  a 
knife-edge  but  a  wedge  which  has  an  angle  at  it’s  apex.  In  the  next  section,  McNamara’s 
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equation  (Eq.3.2)  for  a  knife-edge  will  be  compared  to  the  Balanis’  equation  for  a  wedge. 


3.4  Balanis’  Wedge  Diffraction  Equation 

One  equation  that  calculates  Eiotai  at  the  Rx  by  considering  the  angle  of  the  obstacle 
is  suggested  by  Balanis  [Balanis,  1989]  Balanis  derived  Efotai  has  the  form  of 

Etotai  =  Eo-V-A-e-^^^  (3.9) 

where 

V  =  Diffraction  coefficient 
A  =  Spatial  attenuation 
f3  =  Wave  Number(=  2n/X) 
s  =  Distance  from  the  Obstacle 

V  and  A  are  function  of  the  incidence  angle  (<^’),  the  observation  angle  and  the  wedge 
angle.  Also  V  and  A  are  get  out  of  values  with  different  angle  of  obstacle. 


3.5  Knife-Edge  and  Wedge  Diffraction  Comparisons 

To  verify  that  McNamara’s  equation  is  a  reasonable  approximation  for  the  diffracted 
field  from  wedge  like  obstacles,  the  following  numerical  comparisons  were  performed: 

1.  Compare  Etotai  computed  with  McNamara’s  knife-edge  equation  to  E^otai  computed 
with  Balanis’  wedge  equation  with  the  wedge  angle  set  to  zero. 

•  Estimate  </>,  0’  and  p  for  a  given  obstacle.  Calculate  x  and  y. 

®Balanis  calculates  Etotai  with  a  diffraction  coefficient  and  spatial  attenuation.  In  Appendix  C.  Balanis  ’ 
Wedge  Diffraction  Equation,  whole  equations  are  derived. 
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•  Calculate  Eiotai  with  McNamara’s  knife-edge  equation  (Eq.3.2). 

•  Calculate  Efotai  with  Balanis’  wedge  equation^  with  the  wedge  angle  of  the 
obstacle  set  equal  to  zero. 

•  Compare  the  difference  in  dB  between  \Etotai\^  computed  using  the  McNamara 
and  Balanis  equations. 

2.  Compare  Etotai  computed  with  McNamara’s  knife-edge  equation  to  Etotai  computed 
with  Balanis’  wedge  equation  for  a  variety  of  wedge  angles. 

•  Estimate  and  p  for  a  given  obstacle.  Calculate  x  and  y. 

•  Calculate  Etotai  with  McNamara’s  knife-edge  equation  (Eq.3.2).  In  McNamara’s 
equation,  the  angle  of  the  obstacle  is  ignored  because  the  actual  obstacle  is 
approximated  by  a  knife-edge. 

•  Calculate  Etotai  with  Balanis’  wedge  equation  to  account  for  the  wedge  angle  of 
the  obstacle.  Of  course  Etotai  will  be  different  for  obstacles  with  different  wedge 
angles. 

•  Compare  the  diflference  in  dB  between  \Etotai\'^  computed  using  the  McNamara 
and  Balanis’  equations  for  wedge  angles  ranging  from  to  f  • 

Fig.3.4  shows  examples  of  \Etotai\'^  difference  (dB)  between  the  McNamara  knife-edge  and 
the  Balanis  wedge  computation  {angle  of  obstacle  is  zero).  The  results  for  \Etotai\^  com¬ 
puted  by  the  two  different  methods  are  very  close  unless  (f>  ±  =  nir  where  n  is  an  integer®. 

Etotai  varies  with  the  wedge  angle  of  the  obstacle(wedge)  in  Balanis’  theory.  Fig.3.5 
shows  examples  of  the  difference  (dB)  between  computating  \Etotai\^  with  McNamara’s 
knife-edge  calculation  and  Balanis’  equation  for  a  general  wedge.®  The  results  of  the 
lEtotail"^  comparison  is  very  close  unless  the  wedge  angle  is  large.  Thus  obstacles  which 
roughly  correspond  to  wedges  (e.g.  mountains)  can  be  reasonably  replaced  by  knife-edges 
and  McNamara’s  equation  can  be  used  to  calculate  Etotai  for  many  kinds  of  natural  and 

^See  Appendix  C.  Balanis’  Wedge  Diffraction  Equation. 

^Difference  of  \Etotai\'^  in  decibel  was  calculated  with  a  computer  program(MATLAB).  The  Diffraction 
coefHcient(D)  is  a  function  of  cot{(l>  ±  0’).  If  0  ±  0’  ==  n  ir,  then  cot{(j>  ±  0’)  is  infinite  and  the  computation 
is  invalid. 

®See  Appendix  G.  Balanis  Wedge  Diffraction  Equation. 
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man-made  obstacles.  Therefore,  throughout  this  thesis,  McNamara’s  equation  (Eq.3.2) 
will  be  used  to  estimate  any  additional  path-loss  due  to  diffraction  from  an  obstacle. 


Difference  for  frequency,  phl«1 .25‘pl 


(a),  (f)  =  1.25  TT 


Difference  for  frequency,  phl-1 .15*pi 


— 

X 

— 

-  f-30MHz 

-  f-45MHz 

0 

■  m. 

X 

f-60MHz 

f-IOOMHz 

•  f-200MHz  •  - 

f-300MHz 

• . ■ 

_ _ - 

°  « 

0 

^ ^ ^ ^ ^ ^ ^ ^ - 1 

1  1.2  1.4  1.6  1.0  2  2.2  2.4  2,6  2.8  3 

Observation_ang(radlan)  -  Incidenc6_ang(radian) 


(b).  (j)  =  1.15  TT 


Figure  3.4  The  difference  between  diffraction  computations  with  McNamara’s  knife  equa¬ 
tion  and  Balanis’  general  wedge  equation  for  the  case  of  a  knife-edge  (i.e.  a 
wedge  with  zero  angle),  phi  is  the  observation  angle  {(f>). 


Difference(freq-30MHz),  phi-l  .25*pl 


(a).f  =  30MHz 


Difference(freq-50MHz).  phl-1 .25*pi 


(b).  f  =  50MHz 


Figure  3.5  The  difference  between  the  knife-edge  and  the  wedge  diffracted  fields  for  var¬ 
ious  wedge  angles,  phi  is  the  observation  angle  {(j)). 
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3.6  Summary 

In  this  chapter,  expressions  for  Etotai  based  on  the  geometric  theory  of  diffraction 
or  GTD  were  given.  After  comparing  computations  using  Balanis’  wedge  equation  with 
computations  using  McNamara’s  knife-edge  equation,  one  can  conclude  the  McNamara 
knife-edge  equation  is  a  reasonable  model  (approximation)  for  a  variety  of  wedge-like  ob¬ 
stacles.  Thus,  additional  path-loss  due  to  diffraction  can  be  computed  by  modeling  the 
obstacle  as  a  knife-edge.  The  extra  precision  gained  with  the  Balanis  wedge  equation  is  not 
worth  the  associated  computational  burden.  Hence,  McNamara’s  equation  will  be  used  to 
estimate  additional  propagation  path-loss  due  to  diffraction  throughout  this  thesis. 

Normally,  propagation  path-loss  can  be  calculated  from  Etotai-  Received  signal  power 
(Pr)  has  the  form  of  \Etotai\‘^-  \Etotai\  is  the  magnitude  of  the  total  field  at  the  receiving 
point.  In  the  next  chapter,  propagation  path-loss  prediction  theory  will  be  introduced. 
This  theory  is  the  basis  for  deriving  a  practical  propagation  channel  model  for  electromag¬ 
netic  propagation  over  rough  terrain. 
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IV.  Propagation  Path-Loss 


4.1  Overview 

The  actual  propagation  path-loss  often  exceeds  the  free  space  or  plane  earth  modeled 
path-loss  by  several  orders  of  magnitude.  Furthermore,  it  is  highly  variable  and  fluctuates 
randomly  as  the  Rx  moves  over  irregular  terrain.  Hence,  the  prediction  of  path-loss  is  a 
very  important  step  in  planning  mobile  radio  systems  and  accurate  prediction  methods  are 
needed  to  determine  the  parameters  of  a  radio  system  that  will  provide  efficient  and  reli¬ 
able  service  area  coverage  [Parsons,  1992].  In  this  chapter,  a  path- loss  prediction  method, 
based  on  the  propagation  of  electromagnetic  fields  over  obstacles,  will  be  introduced.  A 
sample  computer  application  will  be  presented  in  section  4.5. 


4.2  Path-Loss  for  One  Obstacle 

An  obstacle  between  the  Tx  and  the  Rx  causes  additional  power  loss  or  path-loss. 
This  diffraction  path-loss  can  take  several  forms.  One,  shadow  diffraction  path-loss  (Lshadow): 
occurs  when  the  Rx  is  in  the  geometric  shadow  region.  Another,  lit  diffraction  path-loss 
(Liit),  occurs  when  the  Rx  is  outside  of  the  shadow  region  and  is  related  to  multi-path 
interference. 


4.2.1  Shadow  Diffraction.  Fig.4.1  shows  a  communication  link  which  has  one 
obstacle  between  the  Tx  and  the  Rx.  If  the  effect  of  the  obstacle  is  ignored  ,  the  received 
signal  power  (P^)  will  be  estimated  using  the  free  space  formula  as  following: 


PtGtGA^ 

(47rd)2 


(4.1) 


However,  if  obstacle  blocks  the  direct  LOS  between  the  Tx  and  the  Rx,  then  a  signal 
will  be  generated  in  the  shadow  region.  Normally,  in  the  shadow  region,  the  magnitude 
of  the  electromagnetic  field  at  the  Rx  is  less  than  the  free  space  field.  Hence,  the  power 
attenuation  (Lghadow)  caused  by  an  obstacle  alters  P^  from  that  predicted  by  Eq.4.1. 
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observation 
angle  {^) 


. / 

i  Incidence; 
;angle((ti’)  ; 


Tx  Edge  Rx 

Figure  4.1  One  obstacle  between  Tx  and  Rx 


Pr  = 


PjGtGrX^ 
(47rd)  ^  L  shadow 


(4.2) 


Lghadow  in  Eq.4.2  is  the  power  loss  in  addition  to  free  space  loss.  Lghadow  can  be  calculated 
from  the  total  electric  field,  Etotaii  at  the  receiving  point  for  the  case  of  the  Rx  in  the 
shadow  region  .  The  power  density  at  the  receiving  point  is 


c  _  PtoialE total 

~  2Z 


(4.3) 


where  Z  is  the  impedance  of  the  medium  and  Etotai  is  given  by  Eq.3.2  which  is  repeated 
here  for  the  reader’s  convenience. 
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Etotai  (2-)  y) 

P'.p-Kkx-7tl4)  1  1 

=  - {l2-C(-7!/)l-jl2-S(-7!/)l} 

(4.4) 

X 

=  p  •  COs[((l>  —  (j)')  —  tt] 

(4.5) 

y 

=  —p  •  sin[{(f)  —  <f>')  —  tt] 

(4.6) 

7 

=  2/\/Ax 

(4.7) 

Eo 

=  Free  space  electric  field 

(4.8) 

Hence, 

5a„  =  ^\{[\-C{-iy)?  +  [\-S{-^y)f){Wm-^)  (4.9) 


or 


5a„  =  S,\{\\-C{-^y)f  +  {\-S{-^y)f){Wm-^)  (4.10) 


where 


5o  = 


EqEI 
2Z  ' 


(4.11) 


where  is  free  space  power  density.  The  unit-less  relative  value  (5rei)  is  defined  as  the 
actual  average  received  power  density  normalized  by  the  free  space  power  density. 

Sr„  =  ^  =  i([|  -  C(-7S)P  +  -  S(-79)P)  (4.12) 

S„  =  (4.13) 

Therefore 

Pr  =  SavAr  (4.14) 

where  Ar  is  Rx  antenna  effective  area  (m^)  and 

.4.  =  ^  (4.15) 

47r 
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(4.16) 


PtGt 

A'Kd? 


Thus 


Pr 


PtGtGrX^ 

(iTTfiy 


1(1^  -  C(-7!/)l’  +  |1  -  S(-79)l"). 


(4.17) 


Or  written  another  way 


P  shadow 


and 


L 


shadow 


1 

m-Ci-iy)?  +  {h-S{-jy)yy 


(4.18) 


(4.19) 


4-2.2  Lit  Diffraction.  Lshadow  occurs  when  an  obstacle  blocks  the  direct  LOS 
between  the  Tx  and  the  Rx.  However,  even  if  the  obstacle  did  not  block  the  direct  LOS, 
the  path-loss  may  still  exceed  the  free  space  case.  This  path-loss  is  called  lit  diffraction 
path-loss  (Liit).  The  term  lit  means  being  excited  by  the  incidence  wave  (field). 

For  this  case,  the  total  field  at  the  Rx  is  composed  of  two  terms.  One  is  the  incidence 
field,  the  other  is  the  diffraction  field.  Thus  the  total  electric  field  at  the  Rx  is  given  by 
Eq.3.2  which  are  repeated  here. 


where 


—  Pincidence  T  Pdif fraction 


7  =  2ly/\x 


(4.20) 

(4.21) 


(4.22) 


The  X  and  y  values  are  given  by  Eq.3.3  and  Eq.3.5  which  are  also  repeated  here  for  conve¬ 
nience. 


X  =  p  •  cos[7r  —  {<f>  —  (f)')] 
y  —  p  ■  sin^K  —  ~  ^')]- 
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It  is  often  convenient  to  express  Lm  in  terms  of  the  first  Fresnel  zone  radius,  Fi,  and  the 
path  clearance,  he-  Referring  to  Fig.4.2,  the  radius  of  first  Fresnel  zone  (Fi)  is  [Flock,  1987], 


(4.23) 


To  avoid  significant  additional  path-loss,  a  clearance  of  about  (0.6)Fi  or  more  is  required. 


EDGE 


Figure  4.2  Subpath  Diffraction 

Defining  the  diffraction  parameter  (ti)  equal  to  y/2  the  total  electric  field  will  be  [Flock. 
1987]; 


Etotat  =  - {[2  +  <^(^)]  -j[^+  5(^;)]} 


(4.24) 


C  — 
— 


(4.25) 


Hence, 


c  — 

^0,7)  - 


^1([1  +  C(»)P  +  |1 
2Z  2  ^2  ‘  4 


([o+C’(u)r  +  [^  +  S(u)]^)(W^m-^) 


(4.26) 
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or 


=  5o^([^  +  C-(r;)f +  [^  +  5W]^)(Wm-2)  (4.27) 


where 


5o  = 


EqEq 

2Z  ■ 


(4.28) 


Hence,  the  unit-less  relative  value  (Srei)  is 

s„,  =  ^  =  +  c(,f  +  [i  +  S(«))")  (4.29) 

and 

=  5o  •  5,e,  (4.30) 


Therefore 


Pr 


^av-^v 


(4.31) 


where  A,,  is  Rx  antenna  effective  area  [rn?)  and 


Ar 

S, 


47r 

PtGt 

Att(P 


(4.32) 

(4.33) 


After  substitution  of  the  above  quantities,  the  received  power  can  be  written  as 


(4.34) 


’■  (AndyLut 

where 

^  M+^w^nF+5(^- 


(4.35) 


(4.36) 
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4-2.3  Multiple  Obstacle  Considerations.  The  analysis  as  described  above 
applies  to  the  case  of  only  one  obstacle  between  the  Tx  and  the  Rx.  However,  in  general 
ground  communication,  there  are  many  obstacles  that  cause  signal  diffraction  and  path- 
loss.  Hence,  a  channel  model  or  method  that  calculates  propagation  path-loss  over  rough 
terrain  must  be  able  to: 

1.  Analyze  terrain  data  to  determine  the  Lit  and  Shadow  diffraction  points. 

2.  Replace  the  diffraction  points  with  a  multiple  knife-edge  sequence. 

3.  Calculate  Lghadow  Lm  for  propagation  over  multiple  knife-edges. 

In  the  next  section,  a  method  for  calculating  path-loss  over  rough  terrain  that  has  multiple 
diffraction  points  (knife-edges)  will  be  introduced. 


4-3  Path-Loss  for  Multiple  Diffraction  Points 

The  diffraction  path-loss  for  multiple  diffraction  points  is  divided  into  two  terms. 
One  is  L shadow:  the  other  is  Lm-  L shadow  may  be  estimated  by  using  a  sequence  of  shadow 
diffraction  knife-edges  (shadow  edge).  However,  Lm  is  not  easily  derived  from  a  sequence  of 
knife-edges.  In  this  section,  a  method  for  computing  Lghadow  foJ^  multiple  diffraction  points 
will  be  introduced  and  a  new  method  that  calculates  Lm  for  rough  terrain  will  be  suggested. 


4-3.1  Shadow  Diffraction.  It  was  suggested  that  a  series  of  obstacles  can  be 
approximated  by  a  sequence  of  knife-edges  as  a  means  of  obtaining  the  path- loss  in  a  simple 
manner.  If  we  find  the  shadow  diffraction  points  that  satisfy  the  following  conditions,  each 
point  can  then  be  approximated  by  a  shadow  edge. 

1.  Each  diffraction  point  must  lie  on  or  above  a  direct  LOS  line  between  the  Tx  and 
the  Rx. 

2.  Each  diffraction  point  lies  along  the  shortest  propagation  path  between  the  previous 
diffraction  point  and  the  next  diffraction  point. 
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3.  The  LOS  between  each  adjacent  diffraction  point  must  be  free  of  obstruction. 

Fig.4.3  shows  how  a  sequence  of  shadow  edges  are  used  to  model  real  terrain.  Note  that 
the  two  diffraction  points  meet  all  three  of  the  conditions  stated  above.  Prom  the  sequence 
of  shadow  edges,  the  total  Lghadow  between  the  Tx  and  the  Rx  can  be  estimated.  Each 
shadow  edge  results  in  a  loss  that  be  estimated  with  Eq.4.19. 


Real  Terrain  between  Tx  and  Rx 


Tx  Rx 


Figure  4.3  Real  terrain  and  sequence  of  shadow  knife-edge 


r  1 

“  |([|  -  ^(-73/)]^ -(- [i  -  5(-72/)]^) 

(4.37) 

X  —  p  •  cos [(</>  — (^')  —  7r] 

(4.38) 

y  =  —p  ■  sin[{(f>  —  4>)  —  tt] 

(4.39) 

7  =  2/\/^ 

(4.40) 

where 


(f)'  :  Incidence  angle  from  the  front  edge  See  Fig.4.1 

(f)  :  Observation  angle  to  the  next  edge  See  Fig.4.1 


The  total  shadow  diffraction  loss,  Lg^adow  (in  dB),  between  the  Tx  and  the  Rx  can  be  cal¬ 
culated  by  simply  summing  each  of  the  individual  Lghadow  (in  dB)  due  to  each  knife-edge 
along  the  path. 
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N 

^shadow  Ei.  N  :  Number  of  shadow  diffraction  edges  (4-41) 

i=l 

4.3.2  Lit  Diffraction.  Eq.3.2  gives  the  total  field  at  the  receiving  point. 

In  the  lit  region,  the  total  field  is  composed  of  an  incidence  field  and  a  diffraction  field.  The 
diffraction  field  due  to  one  subpath  obstacle  can  be  obtained  by  subtracting  the  incidence 
field  from  the  total  field. 


E. 


dif  fraction 


incidence 


(4.42) 


The  total  field  in  the  lit  region  due  to  multiple  subpath  obstacles  is 

N 

Etotal  —  El  -\-  Eincidence  (4-43) 

i=l 

where  N  is  the  number  of  lit  edges  (subpath  obstacles)  and  Ei  is  the  diffraction  field, 
Edif fractions  due  to  the  i-th  lit  edge  (subpath  obstacle).  However,  Eq.4.43  is  not  practical 
for  calculating  Lm  over  rough  terrain  because  of  the  many  possible  path  one  could  con¬ 
sider.  Hence,  a  new  method  that  represents  lit  diffraction  path-loss  over  multiple  subpath 
obstacles  using  path  clearance (h,.)  and  the  first  Fresnel  zone  radius(Fi)  is  simple  to  use. 
For  a  single  subpath  obstacle,  the  lit  diffraction  loss  is  given  by^ 

LiitidB)  =  6- ^-10  for  0  <  h,  <  (0.6)Fi  (4.44) 

El 

Liit{dB)  =  0  otherwise  (4.45) 


For  multiple  subpath  obstacles,  only  the  obstacle  point  that  is  closest  to  the  LOS  (line  of 
sight)  between  Tx  and  Rx  is  used  to  evaluated  and  Fi  in  Eq.4.44.  The  effects  of  all  other 
obstacles  are  ignored.  Fig.  4.4  shows  how  real  terrain  can  be  modeled  for  lit  diffraction 
calculations.  The  parameters  in  Fig.  4.4  are  defined  as  follows: 

•  he'.  Height  difference  between  LOS  (line  of  sight)  and  terrain  point  that  is  closest  to 
the  LOS. 

^See  Appendix  D.  Lut  for  multiple  subpath  edges 
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Terrain  Condition  between  Tx  and  Rx  One  Knife-Edge  between  Tx  and  Rx 


Tx  Rx  Tx  Knife-Edge  Rx 


Figure  4.4  Modeling  real  terrain  with  knife-edges  for  computing  Lm 


•  dt'.  Distance  between  the  Tx  and  the  point  that  is  closest  to  the  LOS. 


•  dr'.  Distance  between  the  diffraction  point  and  the  Rx. 


Radius  of  First  Fresnel  Zone(= 


In  the  next  section,  a  method  for  calculating  path-loss  from  real  terrain  data  will  be 
introduced.  With  this  method,  we  can  create  a  real  terrain  channel  model  that  calculates 
overall  path-loss. 


Real  Terrain  Diflfraction  Model  (RTDMOD) 

In  the  previous  section,  L shadow,  and  Lut  were  derived.  With  these  additional  path- 
losses,  the  received  signal  power,  arriving  over  a  rough  terrain  channel  can  be  estimated 

by 


PtGtGrX^ 

Lshadowl^lit^atmosphere 


(4.46) 


where  Latmosphere  IS  the  additional  path-loss  due  to  the  atmosphere  between  the  Tx  and 
the  Rx.  However,  throughout  this  thesis,  atmospheric  effects  are  neglected.  Hence,  the 
received  signal  power  (P^)  is  given  by 


PtGtGrX^ 

Lshadow^Ut 


(4.47) 
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In  Eq.4.47,  Pr  is  clearly  dependent  on  wavelength  (A)  and  communication  distance  (d) 
which  determine  the  free  space  path-loss  (Lfree)-  Thus  we  can  express  in  the  following 
compact  form: 

Pr  =  (4.48) 


where 


Ltotal 

Lfree 


—  LshadowLutL  free 

_  {Airdy 
A2  ■ 


(4.49) 

(4.50) 


The  total  path-loss  (Ltotai)  between  the  Tx  and  the  Rx  can  be  calculated  from  real  terrain 
data  with  the  following  steps; 


1.  Calculate  Pa/iarfow* 

(a)  Model  the  terrain  as  a  sequence  of  shadow  diffraction  edges  using  the  procedures 
described  in  section  4.3.1. 

(b)  Calculate  Lghadow  for  each  shadow  diffraction  edge  using  Eq.4.37. 

(c)  Calculate  total  Lghadow  by  summing  the  individual  diffraction  losses. 

2.  Calculate  Lm. 

(a)  Find  the  lit  diffraction  edge  location  as  described  in  section  4.3.2. 

•  Consider  only  the  terrain  data  between  the  last  shadow  diffraction  knife- 
edge  and  the  receiver^. 

•  Find  the  point  that  is  closest  to  the  LOS  (line  of  sight)  between  the  last 
shadow  diffraction  knife-edge  and  the  Rx.  The  other  points  will  be  ignored. 

(b)  Calculate  Lm  using  Eq.4.44. 

3.  Calculate  Lfree  using  Eq.4.50. 

^The  lit  diffraction  resulting  from  a  knife-edge  that  is  in  front  of  the  last  shadow  diffraction  knife-edge 
can  be  ignored  because  additional  shadow  diffraction  will  prevent  the  lit  diffracted  field  from  making  a 
significant  contribution  to  the  total  field  at  the  Rx. 
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4.  Compute  Ltotai  (in  dB)  by  summing  the  different  losses.  Ltotai  (dB)  =  L^hadow  (dB) 
+  Liit  (dB)  +  Lfree  (dB) 

Fig. 4. 5  shows  how  real  terrain  data  can  be  replaced  by  an  appropriate  sequence  of  knife- 
edges.  To  simplify  subsequent  discussion,  this  method  of  computing  diffraction  loss  will 
be  called  the  Real  Terrain  Diffraction  Model  (RTDMOD). 


_ 

;  diffraction'. 

i 

Kpdges  y 

1 

I  Find 

(Sequence  of 
i  Sftaefoiv 
\  0//fracf/on 
\Ee/go 


Find 

Lit 

DWractfon 

Edgo 

that  closest 
to  the  LOS 


Tx  Rx 

-  - . ♦  Shadow  diffraction  field 

*  Lit  diffraction  field 


. 

^  /■$hadw\/ 

'!  diffraction', 

,'ru . 't*T 

Kedges  / 

:  diffraction'.  1 

\edpe  J  1 

Figure  4.5  Sequence  of  knife-edges  from  real  terrain 


J^.5  RTDMOD  Computer  Simulation 

The  method  for  calculating  Ltotai  from  real  terrain  data  discussed  in  the  previous 
section  was  implemented  in  a  MATLAB  simulation  (see  Fig.4.6).  In  this  section,  some  of 
the  simulation  results  will  be  discussed. 

The  simulation  depicted  in  Fig.4.6  computes  propagation  path-loss  between  a  specific 
Tx  point  and  each  Rx  point.  Only  the  shortest  communication  path  between  the  Tx  and 
each  Rx  point  was  considered  for  this  simulation. 

A  1km  X  1km  section  of  the  state  of  Texas  as  shown  in  Fig.4.7  was  used  as  the  real 
terrain  data  input.  For  this  example  simulation,  30MHz  frequency  was  used  for  the  carrier 
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(a)JANGl.m  (MATLAB  Function)  (b)PL-total.m  (MATLAB  Function) 


Figure  4.6  Outline  of  simulation 


frequency. 

Lfotai  (dB),  the  total  propagation  path-loss,  is  shown  in  Fig.4.8.  L^hadow,  Lm,  and 


lenBln  Data 


(b)Terrain  Image  Plot 


Figure  4.7  Input  Terrain  data 


Lfree  for  each  receiving  point  are  individually  shown  in  Fig.4.9.  Then  with  available 
for  each  receiving  point  and  knowledge  of  the  transmit  power,  Pj,  and  the  antenna  gain, 
Gt  and  we  can  easily  compute  the  received  power,  Pr,  for  each  point  with  the  following 
dB  version  of  Eq.4.51: 


Pr{dB)  =  Pt{dB)  +  Gt{dB)  +  GridB)  -  LtotaiidB).  (4.51) 
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Figure  4.9  Path-Loss  profile 


4.6  RTDMOD  Validation 

Based  on  Bullington’s  ideas  [Bullington,  1947],  some  techniques  which  predict  the 
power  attenuation  due  to  signal  diffraction  were  developed.  Of  these  techniques,  the 
Epstein-Peterson  (EP)  and  the  Deygout  (DG)  methods  lead  to  the  familiar  GTD  results 
[Pogorzelski,  1982].  Furthermore,  the  path-loss  calculation  model  introduced  in  the  previ- 
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ous  section  is  also  based  on  the  GTD.  In  this  section,  the  validity  of  this  new  model,  the 
Real  Terrain  Diffraction  Model  (RTDMOD),  will  be  confirmed  by  comparison  with  both 
the  EP  and  the  DG  methods. 

4-6.1  Multiple  Edge  Construction  of  EP  and  DG  Method.  In  both  the 
EP  and  DG  methods,  the  diffraction  loss  is  dependent  on  the  distance  parameters  (the  sets 
{s„}  and  {dn}).  These  distance  parameters  are  calculated  using  the  constructions  shown 
in  Fig.4.10.  Notice  that  the  set  {s„}  is  the  same  for  both  the  EP  and  DG  constructions 
but  the  set  {d„}  is  defined  differently  in  each  case. 

Using  these  distance  parameters,  the  propagation  path-loss  can  be  predicted  by  the  EP 


(b).  Deygout  construction 
Figure  4.10  Multiple  Edge  Construction 


and  DG  methods  in  the  following  manner  [Pogorzelski,  1982]: 


•^iotal 


=  ( 
=  ( 


kd\ 


ksiks2 

kdi 


kdo 


ks2ks3 

kdo 


kd: 


ks3ks4 

kds 


ksiks2  ksikss  kssksi 


)  •  [2^/^]-® 


EP  construction 
DG  construction 


(4.52) 

(4.53) 
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where 


k  =  27r/A  (4.54) 

4-6.2  Comparison  Results.  Fig.4.11  compares  EP,  DG,  and  RTDMOD  pre¬ 
dictions  for  a  variety  of  sample  multiple  edge  geometries  over  a  wide  frequency  range.  The 
predictions  of  the  three  GTD  based  methods  differ  by  less  then  5  %. 

However,  both  the  EP  and  the  DG  methods  have  significant  limitations  as  the  num¬ 
ber  of  obstacles  increases.  Even  though,  both  methods  show  good  agreement  with  the 
rigorous  theory  for  a  small  number  of  edges,  the  accuracy  of  both  methods  decrease  as  the 
number  of  obstructions  increase  whereas  the  accuracy  of  RTDMOD  is  not  limited  by  the 
number  of  obstructions. 
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4-7  Summary 

Path-loss  was  derived  using  the  definition  for  the  electric  field  and  power  density. 
Starting  with  Etotai  at  the  receiving  point,  a  path-loss  calculation  method  was  derived. 
RTDMOD  has  the  advantage  of  predicting  the  significant  diffraction  point  locations  as 
well  as  the  total  path-loss  for  a  given  set  of  real  terrain  data.  Unlike  previous  methods, 
this  new  method  can  account  for  an  almost  unlimited  number  of  obstructions.  Although 
RTDMOD  involves  complicated  calculations,  with  the  speed  of  modern  computers,  it  is 
now  possible  to  rigorously  estimate  propagation  path-loss  from  terrain  data.  Nonetheless, 
the  method  described  in  section  4.5  is  computationally  intensive.  Hence,  a  computation¬ 
ally  efficient  stochastic  channel  model  that  predicts  propagation  path-loss  as  a  function 
of  terrain  statistics  (mean,  variance)  would  be  very  useful.  This  is  the  topic  of  the  next 
chapter. 
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V.  Universal  Terrain  Channel  Model  (UTCMOD)  and  Applications 


5.  i  Overview 

Once  the  path-loss  calculation  method  (RTDMOD)  introduced  in  the  previous  chap¬ 
ter  is  implemented  in  a  computer  program,  one  can  readily  use  this  method  to  predict  the 
effect  of  obstacles.  In  this  chapter,  a  stochastic  propagation  channel  model  (Universal  Ter¬ 
rain  Channel  Model)  which  predicts  path-loss  in  decibel  as  a  function  of  terrain  statistics 
(mean,  variance),  signal  wavelength  and  propagation  distance  is  introduced.  In  contrast  to 
RTDMOD  path-loss  predictions,  the  Universal  Terrain  Channel  Model  (UTCMOD)  can  be 
easily  done  on  a  hand-held  calculator.  UTCMOD  will  be  compared  with  the  Hata  model  in 
section  5.6  to  see  how  a  model  based  on  numerical  experiments  (UTCMOD)  compares  to 
a  model  based  on  physical  experiments  (Hata  model).  An  extended  UTCMOD  for  spread 
spectrum  (SS)  signals  will  be  derived  in  section  5.7. 


Propagation  Parameters  and  Assumptions 

In  wireless  ground  communication,  the  propagation  path  between  the  Tx  and  the 
Rx  is  not  trivial.  Unlike  free  space,  rough  terrain  channels  are  random  and  are  not  eas¬ 
ily  analyzed.  Hence,  a  broadly  applicable  rough  terrain  channel  model  which  predicts 
path-loss  effects  of  terrain  between  the  Tx  and  the  Rx  is  urgently  needed  for  improv¬ 
ing  wireless  communication  systems.  A  perfect  channel  model  is  impossible.  However,  if 
sufficient  propagation  path-loss  parameters  can  be  determined,  then  a  broadly  applicable 
channel  model  can  be  created  with  these  parameters.  In  this  thesis,  two  terrain  statistics 
(height  mean  and  variance),  signal  wavelength  and  propagation  distance  are  assumed  to 
be  sufficient  parameters.  I  have  assumed  these  parameters  are  sufficient  for  the  following 
reasons: 

1.  Mean  terrain  height  determines  the  distance  into  shadow  region. 

2.  Variance  of  the  terrain  height  determines  the  number  of  diffraction  points  (knife- 
edges). 
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3.  Signal  wavelength  is  a  fundamental  parameter  for  any  propagation  analysis. 

4.  Propagation  distance  is  also  a  fundamental  propagation  parameter. 

Also  the  following  assumptions  are  used  to  create  a  stochastic  channel  model. 

•  Only  first  order  effects  are  considered  significant.  That  is,  the  effects  of  each  path- loss 
factor  (mean,  variance,  wavelength,  and  distance)  are  independent  of  each  other. 

•  Only  Lghadow  will  be  considered.  Lm  is  less  than  6  dB  over  irregular  terrain. 

However,  the  communication  modes  (e.g.  air-to-ground)  and  the  intervening  terrain  heights 
relative  to  the  height  of  the  Tx  and  the  Rx  should  be  considered  in  predicting  additional 
path-loss  over  rough  terrain.  If  the  terrain  blocks  any  portion  of  the  LOS  path,  then 
diffraction  path-loss  must  be  considered.  Hence,  in  the  next  section,  the  fundamental  def¬ 
initions  of  communication  modes  and  new  definitions  of  communication  links  for  different 
rough  terrain  conditions  will  be  introduced. 


5.3  Communication  Link  Conditions 

Communication  links  are  typically  defined  for  the  four  different  modes  shown  in 
Table  5.1.  In  these  fundamental  definitions  for  communication  links,  the  terrain  conditions 


Communication  Link 

Position  of  the  Tx 

Position  of  the  Rx 

Air-to-Air 

Air 

Air 

Air-to-Ground 

Air 

Ground 

Ground-to-Air 

Ground 

Air 

Ground-to-Ground 

Ground 

Ground 

Table  5.1  Fundamental  Definitions  of  communication  modes 


between  the  Tx  and  the  Rx  are  ignored  and  only  the  positions  of  the  Tx  and  the  Rx  are 
considered.  Hence,  additional  diffraction  path-loss  {L shadow)  might  only  be  expected  in 
the  case  of  Ground-to-Ground  communications.  For  the  other  communication  modes,  one 
might  be  tempted  to  neglect  the  additional  difiraction  path- loss  {L shadow)- 

However,  these  assumptions  are  not  true  in  general.  As  shown  in  Fig. 5.1,  it  can 
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be  assumed  that  some  obstacles  block  the  LOS  between  the  Tx  and  the  Rx  in  Air-to- 
Air  communication  whereas  the  Ground-to-Ground  communication  link  has  a  clean  LOS 
between  the  Tx  and  the  RX.  Thus,  more  inclusive  communication  link  definitions  which 
take  into  account  the  terrain  conditions  between  the  Tx  and  the  Rx  are  required.  In  this 
thesis,  the  definitions  depicted  in  Table  5.2  and  Fig.5.2  for  communication  links  will  be 
used..  Note  that  any  of  the  four  link  conditions  in  Table  5.2  could  occur  for  any  of  the 
four  modes  in  Table  5.1. 


Figure  5.1  The  direct  LOS  and  Obstacle 


Communication 

Link 

LOS(line  of  sight) 
between  Tx  and  Rx 

Height  of  Tx 

Height  of  Rx 

LINK  I 

Lower  than  highest  point 
of  the  Obstacles 

Lower  than  highest  point 

LINK  II 

Blocked  by 
Obstacles 

Higher  than  highest  point 
of  the  Obstacles 

of  the  Obstacles 

LINK  III 

Lower  than  highest  point 
of  the  Obstacles 

Higher  than  highest  point 
of  the  Obstacles 

LINK  IV 

Unblocked  LOS 
between  Tx  and  Rx 

No  Consideration 

No  Consideration 

Table  5.2  Terrain  Condition 

These  new  expanded  communication  link  definitions  (Table  5.2)  provide  the  correct 
perspective  for  analyzing  the  effects  of  rough  terrain  for  any  of  the  communication  modes 
(Table  5.1).  Since  the  additional  diffraction  path-loss  will  differ  for  each  case,  each  link 
condition  (I,  II,  III,  IV)  will  result  in  a  different  stochastic  propagation  model.  One  should 
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Terrain  Condition  For  LINK  I 


Terrain  Condition  For  LINK  II 


(a) 

Terrain  Condition  For  LINK  111 


(b) 


(c) 


Rx  Tx 

(d) 


Figure  5.2  Terrain  Condition  for  LINK  I,  II,  III,  and  IV 


note  that  a  particular  mode  (e.g.  air-to-ground)  could  pass  through  a  variety  of  link 
conditions  (I,  II,  III,  IV)  over  time  as  the  Tx  and/or  the  Rx  move  relative  to  one  another 
and  relative  to  the  intervening  terrain. 

In  the  next  section,  a  stochastic  propagation  model  is  introduced.  The  path-loss 
analytic  prediction  method,  the  Real  Terrain  Diffraction  Model  (RTDMOD),  introduced  in 
chapter  4  will  be  used  to  create  a  stochastic  propagation  channel  model  which  will  be  named 
the  Universal  Terrain  Channel  Model  (UTCMOD).  However,  developing  a  stochastic  model 
(UTCMOD)  from  the  analytic  prediction  method  (RTDMOD)  is  a  little  cumbersome. 
Hence,  the  modeling  method  will  be  explained  first  in  section  5.4  and  then  the  stochastic 
model  for  each  communication  link  condition  (I,  II,  III,  IV)  will  be  introduced  in  section 
5.5. 

5.^  UTCMOD  Development  from  Monte-Carlo  Simulations  of  RTDMOD 

In  chapter  4,  a  computer  program  which  calculates  propagation  path-loss  from  real 
terrain  data  was  introduced.  One  can  use  this  diffraction  based  computer  program  to  create 
a  stochastic  channel  model  that  predicts  path-loss  as  a  simple  function  of  terrain  roughness 
statistics.  The  basic  approach  is  to  select  an  appropriate  parametric  form  for  the  stochas¬ 
tic  model  and  then  estimate  the  coefficients  and  exponents  by  running  numerous  Monte 
Carlo  simulations  using  RTDMOD.  The  first  step  is  to  build  a  Random  Terrain  Genera¬ 
tor  (RTG)  to  create  the  synthetic  random  terrain  required  for  the  Monte  Carlo  simulations. 


5.4-1  Random  Terrain  Generator  (RTG).  The  RTG  must  be  able  to  gener¬ 
ate  synthetic  random  terrain  that  has  some  specified  height  statistics.  The  random  terrain 
data  used  in  the  simulations  were  generated  by  the  random  terrain  generator  (RTG)  shown 
in  Fig.5.3.  The  mean(M<er)  and  variance (Vier)  of  the  terrain  height  can  be  estimated  di¬ 
rectly  from  the  terrain  height(iJfer)  data: 

Mter  =  E{Hter}  mean  of  terrain  height  between  the  Tx  and  the  Rx  (^-l) 

Vter  =  E{Hter  “  Variance  of  terrain  height  between  the  Tx  and  the  Rx(5.2) 


5-5 


MAXLAB  Function 
random  number 
(normal  density) 


output 

Figure  5.3  Random  Terrain  Generator  (RTG) 

where  E{-}  is  the  expectation  operator.  However,  only  the  terrain  height  fluctuations 
relative  to  the  transmitter  height  is  important.  Hence,  the  mean  used  to  generate  the  ter¬ 
rain  is  the  relative  mean  {Mtx)  (i-e.  relative  to  the  height  of  the  Tx)  and  is  calculated  by 
subtracting  the  height  of  the  Tx  (Htx)  from  the  expected  value  of  the  terrain  height  {Mter)- 

Mtx  =  Mter  —  Htx  (5-3) 

Prom  now  on,  the  relative  mean  (Mtx)  and  the  variance  {Vfer)  will  be  used  to  characterize 
the  terrain. 

5.4-2  Path-Loss  Generators.  The  random  terrain  data  with  a  specific  variance 
(Vter)  and  relative  mean  {Mtx)  values  can  be  generated  by  RTG.  As  discussed  in  section 
5.2,  it  will  be  assumed  that  the  mean,  variance,  propagation  distance,  and  wavelength 
(frequency)  can  be  considered  sufficient  path-loss  parameters.  Also  the  effect  of  each  path- 
loss  parameters  are  assumed  to  be  independent  of  each  other.  Hence,  it  is  possible  to 
devise  a  stochastic  model  using  the  following  path-loss  generators. 
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1.  MBPLG  (mean  based  path- loss  generator):  Estimate  path- loss  as  a  function  of  the 
mean  terrain  height. 

2.  MVBPLG  (mean  and  variance  based  path-loss  generator):  Estimate  path-loss  as  a 
function  of  the  mean  terrain  height  and  the  variance  of  terrain  height. 

3.  MVWBPLG  (mean,  variance,  and  wavelength  based  path-loss  generator):  Estimate 
path-loss  as  a  function  of  the  mean  terrain  height,  the  variance  of  terrain  height,  and 
the  signal  wavelength. 

4.  MVWDBPLG  (mean,  variance,  wavelength,  distance  based  path- loss  generator):  Es¬ 
timate  path- loss  as  a  function  of  the  mean  terrain  height,  the  variance  of  terrain 
height,  the  signal  wavelength,  and  the  propagation  distance. 


Figure  5.4  Universal  Terrain  Channel  Model  (UTCMOD)  development  steps 


To  help  understand  how  the  stochastic  model  was  developed,  an  example  for  the  LINK  I 
case  will  be  explained  in  detail. 


5. 4-2.1  MBPLG.  Fig.5.4  shows  the  steps  for  developing  the  stochastic 
model.  The  first  step  is  to  create  the  mean  based  path-loss  equation  using  the  MBPLG. 
Referring  to  Fig.5.5,  the  MBPLG  works  in  the  following  way: 
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MBPLG 


Figure  5.5  Mean  Based  Path-Loss  Generator  (MBPLG) 


1.  Assume  the  very  small  variance  (Vter  =  5m)  so  that  the  variance  based  path- loss  can 
be  neglected. 

2.  Fix  the  propagation  distance  (d=1000m)  and  carrier  frequency  (f=30MHz). 

3.  Select  one  specific  mean  value  {Mtx  =  lOm)  and  generate  random  terrain  data  using 
the  RTG. 

4.  Calculate  Lgkadow  with  the  diffraction  based  computer  program  introduced  in  chapter 

4. 

5.  Repeat  steps  1  ~  4  100  times  using  the  same  mean  (10m),  variance  (5m),  distance 
(1000m),  and  carrier  frequency  (30MHz). 

6.  Calculate  the  average  Lahadow  for  the  specific  mean  value  [Mtx  =  10m). 

7.  Change  the  mean  value  and  repeat  step  1  ~  6  for  different  mean  values.  Create  a 
table  of  the  mean  based  values  for  L shadow  shown  in  Table  5.3. 

8.  Assume  the  mean  based  path-loss  equation  has  the  following  form: 

Lshadow  —  [Mtx]'’^  •  k2  (5.4) 
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Mtx 

J^shadow 

Mean  Based  Model  for  Lghadow 

10m 

shadow^  m— 10 

(10)''!  •  k2 

50m 

(50)'“^  •  k2 

100m 

(100)*i  •  k2 

1000m 

(.1^  shadow)  m— 1000 

(1000)*'!  •  k2 

Table  5.3  Mean  based  Lahadow 


9.  Determine  the  value  of  kl  and  k2  which  minimizes  the  difference  between  Table 
5.3  and  Eq.5.4.  The  variance  of  the  difference  between  Table  5.3  and  Eq.5.4  is  the 
smallest  when  kl=0.5  and  k2=:(0.09).  Even  though  a  bias  exists  as  shown  in  Fig.5.6, 
the  bias  can  be  assumed  to  be  a  distance  (1000m)  based  loss^. 


Result  by  Eq.5.4 

bshadow  (dB) 

=  (Mxxf  ®  0.09 


o  Real  Data 
(Table  5.3) 


Lshadow  (dB) 

-  (Mtx)  0.09 
bias 


Figure  5.6  MBPLG  biased  results 

10.  Rewrite  the  mean  based  path-loss  equation  to  account  for  the  bias. 

Lshadow{dB)  —  i/Mxa;(0.09)  +  bias  (5.5) 


^The  bias  value  which  is  assumed  to  be  a  distance  based  loss 


will  be  explained  in  section  5. 4. 2. 4. 
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5. 4-2. 2  MVBPLG.  The  next  step  is  to  develop  the  mean-variance  based 
path- loss  equation  using  the  MVBPLG.  Referring  to  Fig. 5. 7,  the  MVBPLG  works  in  the 
following  way: 


MVBPLG 

j. - j 


Figure  5.7  Mean- Variance  Based  Path-Loss  Generator  (MVBPLG) 


1.  Select  one  specific  variance  (VJer  =  10m). 

2.  Repeat  step  2  ~  7  of  MRPLG  of  previous  section. 

3.  Change  the  variance  and  repeat  step  1  and  2  for  the  new  variance.  Create  a  table 
of  the  mean-variance  based  values  for  Lghadow  Since  the  MBPLG  already  gave  the 
mean  based  term,  one  can  now  solve  for  the  variance  based  terms  in  Table  5.4. 
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^er 

^shadow 

Mean- Variance  Based  Model 

^shadow 

lOm 

lOm 

(■£a/iadou;)m=10,v=10 

[^10  (0.09)  bias]  -f-  (10)*=^  •  H 

50m 

[VsO  (0.09)  +  bias]  +  (lO)*'^  •  kA 

100m 

i.L  shadow^  m=100^v=l0 

[VlOO  (0.09)  -t-  bias]  -h  (10)*=^  •  ki 

; 

1000m 

(-^a/iadoii;)m=1000,v=100 

[VlOOO  (0.09)  -f  bias]  -b  (10)*^^  •  ki 

100m 

10m 

{.L  shadow^  m=lO,v~XQQ 

[ViO  (0.09)  +  bias]  -f  (100)*®  •  fc4 

50m 

[V50  (0.09)  -h  bias]  -|-  (100)*®  •  H 

100m 

(-^s/iadoiij)m=100,t)=:100 

[VlOO  (0.09)  -f-  bias]  -h  (100)*'^  •  fc4 

lOOOm 

(^a/iarfottj)m=1000,tJ=100 

[VlOOO  (0.09)  +  bias]  -h  (100)*=®  •  ki 

i 

\ 

50,000m 

10m 

(-^a/iad£)iu)m=10,v=50000 

[\/l0  (0.09)  -h  bias]  -t-  (50000)*^  .  fc4 

50m 

(■^a/iarfouj)m=50,v=50000 

[VSO  (0.09)  +  bias]  -h  (50000)'=^  •  fc4 

100,v=50000 

[\/100  (0.09)  -h  bias]  -f-  (50000)*^  •  ki 

1 

1 

mi 

i^shadow)m=1000,v=50000 

[VlOOO  (0.09)  -H  bias]  +  (50000)*®  •  kA 

Table  5.4  Mean- Variance  based  Lg^adow 


4.  Assume  the  {AdB)^  in  Table  5.4  has  the  following  form: 


(AdB),  =  [Fter]*'"  •  fc4 


(5.6) 


5.  Determine  the  value  of  k3  and  k4  which  minimizes  the  difference  between  the  {AdB)^ 
values  in  Table  5.4  and  Eq.5.6.  The  variance  of  the  difference  between  Table  5.4  and 
Eq.5.6  is  the  smallest  when  kl=0.75  and  k2=(0.16). 

6.  Rewrite  the  mean-variance  based  path-loss  equation  in  the  following  way: 

^shadow  {dB)  =  ^/ii^(0.09)  -I-  (v;e^)°-^®(0.16)  +  bias  (5.7) 
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5. 4-2. 3  MVWBPLG.  The  next  step  is  to  develop  the  mean-variance- 
wavelength  based  path-loss  equation  using  the  MVWBPLG.  Referring  to  Fig.5.8,  the 
MVWBPLG  works  in  the  following  way: 

MVWBPLG 


Figure  5.8  MVWBPLG  (mean-variance-wavelength  based  path-loss  generator)  and 
Result 


1.  Select  one  specific  signal  wavelength  (A  =  5m)  . 

2.  Repeat  step  1  ~  3  of  the  MVBPLG  described  in  the  previous  section. 

3.  Change  the  wavelength  and  repeat  the  steps  1  and  2. 

4.  Find  the  functional  relationship  between  Lg^adow  and  the  wavelength.  For  different 
wavelengths,  the  mean  based  value  and  the  variance  based  value  double  as  wavelength 
is  reduced  by  one-tenth. 
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Frequency 

(MHz) 

A 

(m) 

^ter 

(m) 

^1 

J^shadow 

Mean- Variance- Wavelength 
Based  Model  for  Lshadow 

3 

100 

10 

10 

(■^5/iarfou;)A=100,m=10,v=10 

V'10(0.09)(10/100)*® 

-b  (10)°'^®(0.16)(10/100)*=® 

+  bias 

50 

\=100,in—50,v=10 

\/50(0.09)(10/100)''® 

-f-  (10)°-^®(0.16)(10/100)*=® 

+  bias 

100 

55 

55 

* 

■nniii 

55 

55 

■fiTM 

55 

55 

■■ 

50000 

55 

” 

55 

30 

mm 

5? 

55 

55 

55 

60 

55 

55 

55 

55 

150 

2 

55 

55 

55 

55 

300 

1 

10 

10 

(■£4/iarfottj)A=l,m=10,v=10 

v/l0(0.09)(10/l)*5 
+  (10)“-^5(0.16)(10/1)*® 

+  bias 

50 

{L  shadow')  X=l,m=50,v=10 

100 

55 

55 

1 

; 

1000 

55 

55 

100 

55 

55 

55 

; 

* 

I 

50000 

55 

55 

55 

Table  5.5  Mean- Variance- Wavelength  relative  L shadow 


5.  Determine  the  value  of  k5.  Derive  the  mean-variance-wavelength  relative  path-loss 
equation.  (10/A)°'^  accounts  for  the  fact  that  the  loss  values  double  as  the  wavelength 
is  reduced  by  one-tenth. 

^shadow  (dB)  =  ^/]V^(0.09)(10/A)°•^  +  (Vter)°-^®(0.16)(10/A)°'^  +  bias  (5.8) 
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5. 4-2. 4  MVWDBPLG  and  Lahadow  for  LINK  I  case.  The  last  step 
is  to  complete  the  mean-variance-wavelength-distance  based  path-loss  equation  using  the 
MVWDBPLG.  Referring  to  Fig.5.9,  the  MVWDBPLG  works  in  the  following  way: 


MVWDBPLG 


Figure  5.9  MVWDBPLG  (mean- variance- wavelength-distance  based  path-loss  genera¬ 
tor)  and  Result 


1.  Select  one  specific  propagation  distance  (d=5000m). 

2.  Repeat  steps  1  ~  4  of  the  MVWBPLG  described  in  the  previous  section. 

3.  Change  the  distance  and  repeat  the  steps  1  and  2.  Create  a  list  of  the  mean- variance- 
wavelength-distance  based  Lahadow  as  shown  in  Table  5.6. 
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Distance 

(m) 

A 

(m) 

Lter 

(m) 

Mtx 

(m) 

I^shadow 

Mean- Variance- Wavelength-Distance 
Based  Model  for  L.hadow 

500 

100 

10 

10 

iLshadow)d-500,\=W0,m=10,v=10 

\/l0(0.09)(10/100)“-^ 

-F  (10)°-^®(0.16)(10/100)'’-^ 

4-  k6  -1-  (k7)  logio(500/1000) 

; 

1000 

55 

55 

100 

55 

55 

55 

; 

; 

5000 

55 

55 

55 

m 

; 

; 

* 

; 

1 

55 

55 

55 

55 

: 

; 

1 

* 

50000 

>5 

55 

55 

55 

55 

Table  5.6  Mean- Variance- Wavelength-Distance  based  Lghadow 


4.  Assume  the  {bias)d  terms  in  Table  5.6  have  the  following  form: 

{bias)d  =  k6  +  {k7)logio{d/100Q)  (5.9) 

5.  Find  the  value  of  k6  and  k7  which  minimizes  the  difference  between  Table  5.6  and 
Eq.5.9.  The  difference  between  Table  5.6  and  Eq.5.9  is  the  smallest  when  k6=20  and 
k7=13. 

6.  Rewrite  the  mean-variance-wavelength-distance  based  path-loss  equation  in  the  fol¬ 
lowing  way: 

L,hadoxv{dB)  =  y]i^(0.09)(10/A)0-3  -1-  (V(er)'’-^^(0.16)(10/A)“'3  +  20  -f  13log,o{d/1000) 

(5.10) 

Finally,  Eq.5.10  is  the  Universal  Terrain  Channel  Model  (UTCMOD)  for  the  LINK  I  case. 
After  numerous  simulations  and  calculations,  Lg^adow  for  the  LINK  II,  LINK  III  and  LINK 
IV  cases  can  be  derived  in  a  similar  manner.  In  the  next  section,  UTCMODs  for  the  other 
communication  links  (LINK  II,  III,  and  IV)  will  be  introduced. 
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5.4-3  L shadow  for  LINK  II/III  and  LINK  IV  cases.  Using  the  same  method 


as  that  used  for  LINK  I,  the  path-loss  equation  for  the  LINK  II  case  is 


LshadowidB)  =  v^(0.11)(10/A)  -f  (0.16)(10/A)  +  4  +  ISZo^jo (d/1000).  (5.11) 


As  shown  in  Fig.5.10,  LINK  II  is  identical  to  LINK  III  except  for  the  direction  of  propa- 


LINK  II 


LINK  III 


Terrain  Condition 
between 
the  Tx 
and  the  Rx 


Same 
knife-edges 
between 
the  Tx 
and  the  Rx 


Same 

angie  (<)>  -  (t>’) 

for 

each 

knife-edge 


Figure  5.10  The  comparison  of  the  terrain  and  knife-edges  for  LINK  II/III 


gation.  After  making  the  following  observations,  clearly  Lghadow  for  LINK  III  is  the  same 
as  LINK  11.  Both  have  the  same 


1.  The  diffraction  points  for  point-to-point  communication  are  the  same  regardless  of 
the  propagation  direction. 

2.  As  seen  in  Fig.5.10,  each  diffraction  point  (knife-edge)  has  the  same  angle  difference 
((/>  -  (j)')  for  both  link  conditions  (LINK  II  and  LINK  III). 

3.  As  explained  in  chapter  3,  diffraction  path-loss  can  be  predicted  from  the  distance 
between  the  knife-edges  (p)  and  the  angle  difference  {(j)  -  0’).  Hence,  L shadow  for 
LINK  III  can  be  estimated  with  same  equation  for  LINK  11. 
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Therefore,  Lghadow  for  LINK  III  is 

L,hadow{dB)  =  x/]i^(0.11)(10/A)  +  (Vier)‘’-'"(0.16)(10/A)  +  4  +  13/0510  WlOOO)  (5.12) 
which  is  identical  to  Eq.5.11. 

If  the  communication  link  between  the  Tx  and  the  Rx  is  the  LINK  IV  condition,  then 
the  LOS  between  the  Tx  and  the  Rx  is  unobstructed  within  at  least  0.6  Fi  and  diffraction 
is  negligent.  Hence,  Lshadow  for  LINK  IV  is  zero. 

Lshadow{dB)  =  0  (5.13) 

5.5  UTCMOD  Summary  and  Sample  Applications 

In  summary,  UTCMOD  for  predicting  the  additional  path-loss  due  to  dififraction  can 
be  written  in  the  following  form. 

LshadoUdB)  =  M{dB)  +  V{dB)  +  R{dB)  (5.14) 

where 

M (dB)  :  mean  based  path- loss. 

V{dB)  :  variance  based  path-loss. 

R(dB)  :  distance  based  path- loss. 

and  the  wavelength  based  path- loss  is  included  in  M  (dB)  and  V  (dB).  The  specific  ex¬ 
pression  for  M,  V,  and  R  for  each  link  case  are  summarized  in  Table  5.7. 
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Lshadou,{m  =  M(dB)  +  V(dB)  +  R(dB) 


Communication 

Link 

M(dB) 

V(dB) 

R(dB) 

LINK  I 

.  (^)0.3  .(0  09) 

i\VTer\r^  ■  ar  -(o.ie) 

20  -1-13  •  logio 

LINK  II/III 

(I^Terl)”-^^  •  (f )  -(0.16) 

4  +13  •  logio 

LINK  IV 

0 

0 

0 

Table  5.7  UTCMOD  parameters 


5.5.1  Frequency  Based  Form  for  UTCMOD.  Table  5.7  can  be  revised  to 
have  expressions  in  terms  of  frequency  rather  than  wavelength  by  noting  that 
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(5.15) 


and 

d 

13  •  =  13  ■  i^ogiod  -  3]  =  13  •  logiod  -  39.  (5.16) 

The  results  of  the  conversion  from  wavelength  to  frequency  are  shown  in  Table  5.8. 


5.5.2  Sample  Application  Using  UTCMOD.  In  the  previous  section,  the 
Universal  Terrain  Channel  Model  (UTCMOD)  for  various  communication  links  was  in¬ 
troduced.  In  this  section,  a  MATLAB  subroutine,  JANG2.m,  that  calculates  additional 
diffraction  based  path-loss  using  the  UTCMOD  will  be  discussed.  The  simulation  depicted 
in  Fig.5.11  computes  Lghadow  between  the  Tx  (origin  point)  and  every  possible  receiving 
point.  L shadow  was  calculated  using  the  following  steps: 


1.  Compute  the  Mtx  and  Vter  between  the  Tx  and  the  Rx  from  the  real  terrain  data. 
Fig.5.12  shows  the  input  terrain  data. 
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Lshado,.{.<^^)  =  M(dB)  +  V(dB)  +  R(dB) 


Condition 

M(dB) 

V(dB) 

R(dB) 

I 

(|^Ter|)°-""  •  (9.14  XlO-4)  • 

-19  -\-\3’logiQd 

II/III 

(ll^Terl)'’-'®  •  (5.33  XlO-9)  .  f 

-35  -fl3.fo5'iod 

IV 

0 

0 

0 

Table  5.8  UTCMOD  parameterized  on  frequency 


2.  Determine  the  Link  condition  between  the  Tx  and  the  RX  and  use  the  appropriate 
equation  from  Table  5.8  to  calculate  Lshadow  The  Link  conditions  were  defined  in 
section  5.3  (Table  5.2). 

3.  Plot  the  results.  Fig.5.13  shows  the  resulting  Lg^adow  predictions  for  the  terrain  data 
shown  in  Fig.5.12. 

5.6  UTCMOD  Validation 

The  prediction  of  path-loss  is  a  very  important  step  in  planning  an  efficient  ground 
wireless  radio  system.  Hence,  a  number  of  prediction  methods  (models)  have  been  de¬ 
veloped.  These  are  exemplified  by  the  Egli  model  \Egli,  1957],  the  Longley-Rice  model 
[Longley  and  Rice,  1968],  the  Okumura  method  [Okumura  et  al.,  1968],  and  the  Hata 
model  [Hata,  1980],  which  are  all  statistical  propagation  models.  The  Hata  model,  which 
is  based  on  the  Okumura  method,  is  a  widely  used  propagation  model  [Delisle  et  al.,  1985] 
[Peterson  et  al.,  1995].  Hence,  in  this  section,  path-loss  predictions  of  my  stochastic  model 
(UTCMOD)  will  be  compared  with  the  predictions  of  the  Hata  model  as  a  preliminary 
validation  of  UTCMOD. 
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Hetght(m) 


L... ................ ............ .T.V........... - * 

(a).  Jang2.m  (MATLAB  Function)  (b).  PL-sto.m  (MATLAB  Function) 


Figure  5.11  Ditfraction  path-loss  compution  using  UTCMOD 


Figure  5.12  Sample  input  Terrain  Data 


Dlstance(m) 


0  0 


Dlstanc6(m) 


4000  6000 

Distance(m) 


(a)  (b) 

Figure  5.13  Lghadow  for  each  receiving  point  (f  =  30MHz) 


5.6.1  Hata  Model.  Hata  presented  the  urban  area  propagation  loss  as  a  stan¬ 
dard  formula  and  supplied  correction  equations  for  other  situations.  The  standard  formula 
for  median  path-loss  in  large  urban  areas  is  given  by 


L{dB)  =  69.55-f26.16Zo5(/e)-13.82/o5(/i<e)-a(/ire)-(44.9-6.55/op(/i(e))+/o5(d)  (5.17) 


where  the  correction  parameter,  a,  is  used  to  reduce  the  mean  error.  For  example, 

(^(Ke)  =  8.29[Zo5'(1.54/ire)]^  “  l.l(dS)  for  large  city,  fc  <  300MHz 

a{hre)  =  3.2[log{11.75hre)]^  —  4:.97{dB)  for  large  city,  fc>300MHz 

In  Eq.5.17,  Hata  model  predicts  total  propagation  path-loss  {Ltotai  =  f'/ree  +  Bskadow) 
as  a  function  of  carrier  frequency  (fc),  Tx  antenna  height  (/ije),  Rx  antenna  height  (hre), 
and  propagation  distance  (d).  The  Hata  model  is  based  on  the  basic  conditions  of  the 
Okumura  method  which  can  be  summarized  as  follows  [Delisle  et  al.,  1985]: 

1.  Propagation  losses  are  computed  between  the  antennas. 

2.  Smooth  terrain  is  assumed  and  the  additional  losses  due  to  each  obstacle  are  not 


taken  into  account. 


3.  A  basic  formulation  for  propagation  losses  is  established  for  an  urban  environment 
and  then  a  correction  factor  must  be  introduced  to  adjust  for  other  cases. 

In  the  Hata  model,  propagation  path-loss  is  very  sensitive  to  propagation  distance  but  the 
terrain  roughness  is  neglected.  Hence,  Hata  model  is  well  suited  for  plane  urban  area  but 
is  insufficient  for  irregular  terrain. 


5.6.2  UTCMOD.  UTCMOD  predicts  additional  path-loss  based  on  signal 
diffraction.  The  model  calculates  path-loss  from  terrain  roughness  statistics  (mean  and 
variance)  and  fundamental  parameters  (distance  and  carrier  frequency).  As  an  example. 
Table  5.9  shows  the  path-loss  parameters  for  LINK  I,  H,  III,  IV  conditions  using  UTCMOD. 

^total  (dB)  =  Lf,.,,idB)  +  L 

shadow  (dB)  (5.18) 

where 

A'jr  ,  (1 

LfreeidB)  =  20logio[—^]  (5.19) 

and 

I^shadow  (dB)  =  M{dB)  +  V{dB)  +  R{dB)  (5.20) 


Condition 

M(dB) 

V(dB) 

I 

ViMr.l  •  (5.14  xlO-^)  .  /o-s 

(I^Terl)”  ""  •  (9.14  XlO-4)  .  f 

-19  -M3-/o5riod 

II/III 

(I^Terl)"-'®  •  (5.33  XlO-®)  •  f 

-35  -|-13-/o5iod 

IV 

0 

0 

0 

Table  5.9  UTCMOD  parametric  terms  for  LINK  I,  II,  III,  IV  conditions 


5-22 


Mainly,  terrain  roughness  (irregularity)  determines  terrain  statistics  which  in  turn 
determine  the  additional  path-loss.  Propagation  distance  contributes  very  little  to  the 
additional  path-loss  because  over  the  rough  terrain,  diffraction  is  the  primary  phenomena. 
However,  distance  is  accounted  for  in  the  free  space  path-loss  term.  Finally,  UTCMOD 
covers  all  communication  link  conditions  and  is  suitable  for  rough  mountain  areas. 

5.6.3  Comparison  of  UTCMOD  and  Hata  model  .  The  Hata  model  and 
UTCMOD  are  based  on  different  terrain  conditions.  The  Hata  model  applies  primarily  to 
urban  areas  whereas  UTCMOD  applies  to  rough  mountainous  terrain.  Also,  the  models  use 
different  parameters  to  predict  path-loss.  Hence,  it  is  of  interest  to  evaluate  how  UTCMOD 
compares  with  the  Hata  model  for  a  specific  scenario.  The  details  of  the  comparison  are 
given  below: 

1.  Hata  Model 

•  Tx  Antenna  height  (hte)  is  the  base  radio  station  antenna  height  in  meters.  Two 
different  hte  (Sm  and  30m)  were  used. 

•  Rx  Antenna  height  (/ire)  is  the  mobile  antenna  height  in  meters.  An  hre  of  Im 
was  used. 

•  Terrain  is  assumed  to  be  a  10m  average  building  height  in  an  urban  area. 

2.  Universal  Terrain  Channel  Model  (UTCMOD) 

•  Tx  relative  mean  (Mtx)  is  the  average  height  of  the  terrain  relative  to  the  height 
of  Tx  point. 

•  Terrain  Variance  (VVer)  is  the  variance  of  the  terrain  height. 

3.  Comparison 

•  For  hte  =5m,  an  average  building  height  of  10m,and  hre  is  Im,  the  points  of 
the  Tx  and  the  Rx;  are  lower  than  terrain  height.  This  represents  the  LINK  I 
condition  in  UTCMOD.  Hence,  in  this  case,  the  Hata  model  is  comparable  to 
the  LINK  I  case  of  UTCMOD. 
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•  For  hte  =30m,  the  point  of  the  Tx  higher  and  the  Rx  is  lower  than  the  terrain 
height  which  represents  the  LINK  II  condition  in  UTCMOD.  Hence,  in  this 
case,  the  Hata  model  is  comparable  to  the  LINK  II  case  of  UTCMOD. 

The  comparison  parameters  are  summarized  in  Table  5.10.  Path-loss  difference  in  deci- 


hte  (m) 

Hata  Model 

UTCMOD 

5 

hre  =  Im 

Average  building  height  :  10m 

Mtx  =  (Average  terrain  height  -  hte)  =  (10m-5m)  =  5m 
Om  <  Vrer  <  100m 

30 

Mtx  =  (Average  terrain  height  -  hte)  =  (10m-30m)  =  -20m 
Om  <  Vxer  <  100m 

Table  5.10  Pertinent  parameters  for  comparing  the  Hata  model  and  UTCMOD 


bels  as  a  function  of  frequency  is  shown  in  Fig.5.14.  UTCMOD  predicts  a  slightly  higher 
path-loss  than  the  Hata  model  because  the  Hata  model  assumed  a  smooth  terrain  and 
the  additional  losses  due  to  each  obstacle  are  not  taken  into  account.  The  models  are  in 
good  agreement  (within  IdB/km)  over  both  a  broad  carrier  frequency  range  and  a  large 
terrain  height  variance  range.  With  UTCMOD,  we  can  also  predict  how  SS,  other  wide 
bandwidth,  and  variable  frequency  communication  schemes,  propagate  over  rough  terrain. 
In  the  next  section,  a  SS  propagation  path-loss  channel  model  will  be  introduced. 


5.7  Spread  Spectrum  Communication  Channel  Model 

Although  SS  methods  are  often  effective  in  reducing  multi-path  interference  problems, 
path-loss  of  SS  in  extraordinarily  complex  propagation  paths  such  as  those  characterized 
by  extremely  rugged  terrain  is  not  well  understood[  Woener  et  al,  1995].  Rough  terrain 
has  many  different  communication  obstacles  which  increase  the  complexity  of  the  direct 
path.  Since  prediction  of  path-loss  is  an  important  step  towards  improving  SS  methods  in 
ground  communication,  a  propagation  channel  model  for  SS  will  now  be  introduced. 


5.7.1  Spread  Spectrum  Channel  Bandwidth.  As  described  in  section  2.3, 
the  transmitted  bandwidth  of  a  SS  signal  is  determined  by  a  pseudo-noise  (PN)  code. 
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Figure  5.14 


Difference  between  the  Hata  and  UTC  models’  path-loss  prediction. 
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In  DSSS  system,  the  bandwidth  is  a  direct  function  of  the  chip  rate.  In  FHSS  system, 
bandwidth  is  determined  by  the  number  of  discrete  hopping  frequencies.  A  SS  system  is 
a  system  in  which  the  transmitted  signal  is  spread  over  a  frequency  band  that  is  much 
wider  than  the  minimum  bandwidth  required  to  send  the  signal.  Hence,  the  propagation 
frequency  (fp)  parameter  for  diffraction  path-loss  of  a  SS  signal  is  different  from  a  non-SS 
signal  (narrow  band  signal)  carrier  frequency  parameter. 

In  the  DSSS  technique,  first  the  PN  code  is  modulated  onto  the  information  signal 
using  one  of  several  modulation  techniques  (eg.  BPSK,  QPSK,  etc  ).  Then,  a  doubly 
balanced  mixer  is  used  to  multiply  the  RF  carrier  and  PN  modulated  information  signal. 
This  process  causes  the  RF  signal  to  be  replaced  with  a  very  wide  bandwidth.  That  means 
there  is  a  broad  band  of  propagating  frequencies  in  addition  to  the  RF  carrier  frequency 
and  the  range  is(iS'ee  the  section  2.3). 


MHz)  <  fp{Hz)<MHz) 

(5.21) 

MHz)  =  MHz)-MiHz) 

(5.22) 

MHz)  =  MHz)  +  M{Hz) 

(5.23) 

where  is  carrier  frequency,  Rc  is  the  chip  rate,  /j  is  the  lowest  frequency,  is  the  highest 
frequency. 

However,  in  FHSS,  the  incoming  digital  stream  is  shifted  in  frequency  by  an  amount 
determined  by  a  code  that  spreads  the  signal  power  over  a  wide  bandwidth.  The  FHSS 
transmitter  is  a  pseudo-noise  PN  code  controlled  frequency  synthesizer.  The  instantaneous 
frequency  output  of  the  transmitter  jumps  from  one  value  to  another  based  on  the  pseudo¬ 
random  input  from  the  code  generator.  Varying  the  instantaneous  frequency  results  in 
an  output  spectrum  that  is  effectively  spread  over  the  range  of  frequencies  generated.  In 
FHSS,  a  number  of  discrete  frequencies  will  be  used  to  the  propagated  signal.  Hence,  the 
range  of  the  propagation  frequency  {fp)  is 

fi{Hz)  <  fp{Hz)  <  MHz)  (5.24) 
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where 

MHz)  -  MHz)  (5.25) 

and 

MHz)  =  [fo  +  k-Af]{Hz)  (5.26) 

where  k  is  the  number  of  discrete  frequencies  that  is  determined  by  the  SS  code  and  Af  is 
the  hop  frequency  spacing. 

The  propagation  frequency  (/p)  of  a  SS  signal  is  different  from  that  of  the  narrow 
band  signal  and  has  the  range  (/;  <  fp  <  fh)-  Therefore,  UTCMOD  (Table  5.8)  should  be 
extended  to  be  a  function  of  /;  and  fh.  to  predict  the  path-loss  of  a  SS  signal.  In  the  next 
section,  a  channel  model  which  estimates  additional  diffraction  path-loss  of  a  SS  signal  will 
now  be  presented. 


5.7.2  A  Modified  UTCMOD  for  Spread  Spectrum  Signals.  Additional 
propagation  path- loss  {L shadow)  can  be  estimated  with  the  relative  path-loss  parameters  ( 
a,  M,  V,  and  R  )  which  were  given  in  Table  5.8.  As  a  example,  Lghadow  for  the  case  of 
LINK  I  is; 

Lghadow  =  y|Mrx|-(5.14xlO-^)./«-3  +  (|Urer|)°-'®-(9.14xlO-'‘)./°-3  +  13-Zo5iod-19  (5.27) 

If  we  assume,  in  SS  communication,  the  probability  of  occurance  for  each  frequency  between 
fi  and  fh  is  the  same,  then  the  additional  path-loss  for  SS  {LSSghadow)  in  decibels  can  be 
approximated  by  the  mean  value  of  Lghadow  in  decibel  over  the  bandwidth  of  the  SS  signal. 


LS  S  ghadowi.^H) 

ll;'[Lghadow{dB)]df 

[fh  ~  fl] 

=  J\MtA  ■  (5.14  X  10-^)  ■  y:  -H  i\VTer\)°-^^  '  (9.14  X  10-^)  ■  +  13  •  log,,d  -  19 

[fh  -  fl\  [fh  -  fl\ 


=  •  (5.14  X  10-^)  •  +  {[VTer[)°  '^^  '  (9.14  X  10-^)  •  //tu.  +  13  •  log,od  -  19 


[fi-^  -  in 

(i.3)[A-/,] 
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Therefore  LSS shadow  for  the  case  of  LINK  I: 


LSSshadow  -  V|MT.|-(3.95xlO-^)-i% - ^+(|K^,,|)°'^^(7.03xl0-^) 

^  [Jh  -  Jl\ 


lfl.3  _  fl.3] 

-4^  ^ +13.logiod-l9 


[fh  -  fi] 


(5.28) 


Using  the  same  method,  additional  propagation  path-loss  for  different  communication  links 
can  be  derived.  The  results  are  given  in  Table  5.11. 

Clearly,  UTCMOD  for  SS  signals  is  different  from  that  for  narrowband  signals  because 


LSSshadow  =  M  (dB)  -k  V  (dB)  -f  R  (dB) 


Link 

condition 

M(dB) 

V(dB) 

R(dB) 

LINK  I 

(|^^^^!)0.r5  .  r_g^  .  [^1.3  .  ^i.3] 

-19  +  ISlogiQ  d 

LINK  II  /  III 

(|p.^^^|)0.r5  .  .  [/2  .  y2] 

-35  H"  TiSlogiQ  d 

IV 

0 

0 

0 

Table  5.11  Shadow  Diffraction  Path-Loss  for  SS 


SS  signals  have  a  range  of  propagating  frequencies.  Also  the  frequency  range  of  DSSS  is 
different  from  FHSS.  However,  the  diffraction  path-loss  effected  by  the  terrain  does  not 
vary  much  over  typical  SS  bandwidths.  Fig.5.15  shows  the  additional  path-loss  results  for 
different  propagating  frequencies  over  the  same  terrain  {See  Fig.  5. 12  for  the  input  terrain 
data  and  Table  5.12  for  the  BW  of  each  of  the  three  cases).  The  results  for  all  three  cases 


Figure 

Propagation  Frequency 

Bandwidth 

Fig.5.15.  (a) 

fp  =  30MHz 

Narrow  Band  Signal 

Fig.5.15.  (b) 

28MHz  <  fp<  32MHz 

BW=4MHz 

Fig.5.15.  (c) 

30MHZ  <  fp<  80MHz 

BW=50MHz 

Table  5.12  Propagation  Frequency  and  Results 


are  very  similar  with  only  a  change  in  the  path-loss  magnitude  scale.  The  purpose  of  this 
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comparison  is  show  the  effect  of  the  terrain  on  SS  signals  and  show  that  knowledge  of 
the  terrain  is  important  toward  understanding  SS  propagation  over  rough  terrain.  Even 
though  SS  methods  are  effective  in  reducing  multi-path  interference  problems,  SS  methods 
are  susceptible  to  signal  diffraction  path-loss.  However,  the  diffraction  path-loss  can  be 
predicted  by  UTCMOD  for  SS  (Table  5.11)  and  prediction  of  such  path-losses  will  be  an 
important  step  towards  improving  SS  methods  in  ground  communication. 


Figure  5.15  Lg^adow  for  each  receiving  point. 


5.8  Summary 

Universal  Terrain  Channel  Model  (UTCMOD),  a  reasonable  propagation  model  for 
estimating  the  mean  path-loss  of  signal  propagation  over  rough  mountainous  terrain,  was 
derived  from  analyzing  the  terrain  effects  on  signal  propagation  and  is  valid  over  various 
frequency  ranges.  UTCMOD  is  applicable  to  many  different  kinds  of  terrain  conditions 
whereas  the  Hata  model  is  most  relevant  for  terrain  conditions  characteristic  of  large  cities 
or  medium-small  cities.  By  considering  communication  link  conditions,  UTCMOD  pro¬ 
vides  the  correct  perspective  for  analyzing  propagation  effects.  In  section  5.6,  UTCMOD 
is  compared  with  the  Hata  model.  The  models  are  in  good  agreement  (within  1  dB/km) 
over  both  a  broad  carrier  frequency  range  and  a  large  terrain  height  variance  range.  In 
section  5.7,  UTCMOD  for  SS  signals  was  derived  because  prediction  of  path-losses  an  im¬ 
portant  step  towards  improving  SS  methods  in  ground  communication. 

In  the  next  chapter,  an  application  of  the  Universal  Terrain  Channel  Model  (UTC¬ 
MOD)  to  the  analysis  of  tactical  communication  will  be  introduced. 
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VI.  An  Application  of  UTCMOD  to  the  Analysis  of  Tactical 

Communication 

6. 1  Overview 

A  systems  engineer  is  interested  in  knowing  the  area  coverage  associated  with  a 
specific  system  in  a  particular  environment.  In  military  tactical  communication,  an  ad¬ 
ditional  objective  is  to  reduce,  as  much  as  possible,  the  range  at  which  the  signal  can  be 
intercepted  by  unfriendly  receivers  while  at  the  same  time  maintaining,  or  even  improving, 
the  communication  range  of  the  intended  receiver.  This  is  referred  to  as  Low  Probability 
of  Interception  (LPI)  communications. 

Existing  LPI  theory  ignores  terrain  conditions  while  other  performance  metrics  of 
communication  systems  are  emphasized  to  satisfy  the  objective  of  LPI  communication. 
However,  as  discussed  in  the  previous  chapter,  terrain  conditions  are  an  important  param¬ 
eter  which  determines  the  received  power,  in  a  wireless  communication  system.  Hence, 
terrain  conditions  should  be  considered  in  the  analysis  of  LPI  communication  systems. 

In  this  chapter,  the  terrain  quality  factor  (Qter)  which  can  be  computed  with  UTC¬ 
MOD  and  sample  applications  will  be  introduced.  First,  basic  LPI  theory  and  the  current 
LPI  quality  factor,  Qlpu  will  be  explained  in  section  6.2  and  section  6.3.  Then  the  terrain 
quality  factor  (Qter)  which  accounts  for  terrain  effects  will  be  introduced  in  section  6.4. 
Finally,  the  tactical  quality  factor  (Qtac)  which  can  be  used  to  estimate  the  probability 
of  interception  and  other  tactical  applications  of  this  new  quality  factor  will  be  introduced 
in  section  6.5. 

6.2  Basic  LPI  Theory 

Basic  LPI  communication  theory  was  developed  to  analyze  typical  military  wireless 
communication  scenarios  and  communication  links  to  assess  their  vulnerability  to  intercep¬ 
tion.  An  overview  of  this  relatively  new  theory  is  now  presented  to  provide  the  necessary 
background  for  understanding  the  changes  suggested  in  subsequent  sections. 
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6.2.1  Typical  Military  Wireless  Communication  Scenario  .  Fig.6.1  shows 

a  typical  military  wireless  communication  scenario.  While  the  communication  Rx^is  the 
intended  receiver,  the  interception  Rx^  is  the  unintended  receiver  which  attempt  to  detect 
and  exploit  the  transmitted  signal  to  intercept  the  transmitted  information. 

The  objective  of  any  LPI  communication  system  is  to  relay  information  from  the 


Figure  6.1  The  typical  military  wireless  communication  scenario 

Tx  to  the  Rx  while  minimizing  the  ability  of  an  unauthorized  listener  to  intercept,  classify, 
or  otherwise  exploit  the  transmitted  signal.  Communication  systems  can  use  a  variety  of 
techniques  for  reducing  the  probability  of  interception  such  as  steerable  high  gain  antennas, 
adaptive  transmitter  power  control,  and  waveforms  with  large  time-bandwidth  products 
and  noise-like  spectra  to  name  a  few.  Likewise,  the  interceptor  has  similar  technologies, 
such  as  directional,  low  sidelobe  antennas  and  adaptive  filtering  [Mills,  1994]. 


6.2.2  Communication  Link  Analysis.  If 

the  received  signal  power  (P^)  is 


PtGtGrX^ 

(47rd)^L 


terrain  conditions  are  ignored,  then 


(6.1) 


where  L  is  the  sum  of  the  all  other  losses  except  L shadow 

In  LPI  communication,  this  simple  link  equation  can  be  used  to  compare  the  received 


^Communication  Rx  means  intended  (or  friendly)  Rx. 
^Interception  Rx  means  unintended  (or  unfriendly)  Rx. 


6-2 


power  at  the  intended  Rx  (Sc)  to  the  received  power  at  the  interception  Rx  (5/).  Thus, 
referring  to  Fig.6.1,  the  received  power  at  the  intended  Rx  (Sq) 


PiGtcGct^^ 

{AnRcYLc 


(6.2) 


where 

•  Sc  is  the  received  signal  power  at  the  intended  Rx 

•  Gtc  is  the  Tx  antenna  gain  in  the  direction  of  the  intended  Rx 

•  Got  is  the  intended  Rx  antenna  gain  in  the  direction  of  the  Tx 

•  Rc  is  the  distance  between  the  Tx  and  the  intended  Rx 

•  Lc  is  the  atmospheric  loss  factor. 

Consequently  the  received  signal  power  to  noise  PSD  can  be  expressed  in  the  terms  of  the 


link  parameters: 


Sc  _  PtGrcGcT ,  A  ,2 

LcNsc  ‘i'^Rc 


N. 


sc 


(6.3) 


Solving  Eq.6.3  for  the  maximum  communication  range  (Rc)  yields 


PiGtcGct  f  A 
Rc  =  \l-  r  Ar - iT-r- 


LcNsc  47r  Sc/Nsc 


(6.4) 


The  maximum  interception  range  (Ri)  can  be  estimated  by  the  same  type  link  pa¬ 
rameters.  Thus  referring  to  Fig.6.1,  the  received  power  at  the  interception  Rx  (S/)  is 


Si  = 


PiGtiGit^ 

{A-nRiTLi 


(6.5) 


where 


•  Gti  is  the  Tx  antenna  gain  in  the  direction  of  the  unintended  Rx 

•  Gjx  is  the  unintended  Rx  antenna  gain  in  the  direction  of  the  Tx 

•  Rj  is  the  distance  between  the  Tx  and  the  intercept  Rx 

•  Li  is  the  atmospheric  loss  factor  between  the  Tx  and  the  intercept  Rx. 
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Analogous  to  Eq.6.3,  the  received  signal  power  to  noise  PSD  ratio  can  be  expressed  in  the 
terms  of  the  link  parameters: 


Si  _  PiGtiGit  ,  ^  >.2 
Nsi  LjNsi  A'kRj 


Solving  Eq.6.6  for  the  maximum  interception  range  (i?/)  yields 


Ri 


PtGxiGiT,  A  1 
LjNsi  47r  Sj/Nsi 


(6.6) 


(6.7) 


6.2.3  Possible  Communication  and  Interception  Areas.  In  the  previous 
section,  expressions  for  the  maximum  communication  range  (Re)  and  the  maximum  inter¬ 
ception  range  (i?/)  were  derived.  Rc  and  Rj  indicate  a  geometric  range  from  the  Tx  and 
thus  determine  the  possible  communication  area  {Ac)  and  the  possible  interception  area 
(A/)  as  shown  in  Fig. 6. 2.  Usually,  the  interception  Rx  requires  a  lower  signal-to- noise 


(a) 


(b) 


Figure  6.2  The  relationship  between  the  geometrical  range  and  the  possible  geometrical 
area.  (a).  Rc  and  the  possible  communication  area.  (b).  Ri  and  the  the 
possible  interception  area. 


ratio  (SNR)  to  detect  the  presence  of  a  signal.  Hence,  without  employing  LPI  techniques, 
Ri  is  usually  longer  than  Rc  and  hence  A/  is  larger  than  Ac.  That  means  the  interception 
Rx  can  detect  the  signal  outside  of  Ac-  However,  when  LPI  techniques  are  employed, 
the  interception  Rx  needs  to  move  closer  to  detect  the  signal.  Some  of  these  techniques 
are  high  gain  directional  antennas,  adaptive  interference  suppression  filters,  and  custom 
wave  design.  The  main  objective  of  these  techniques  and  LPI  communication  is  to  reduce 
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Ri  while  simultaneously  maintaining,  or  even  improving  (i.e.  lengthening),  Rc-  In  other 
words,  the  objective  of  the  LPI  communication  is  increase  the  ratio  of  Rq  to  Ri  {Rq  / 
Ri). 

In  the  initial  LPI  formulation  [Gutman  and  Prescott,  1989],  this  ratio  was  the  pri¬ 
mary  metric  for  assessing  LPI  performance.  Thus  20log{^)  became  known  as  the  LPI 
quality  factor  or  Qlpi-  First,  the  existing  Qlpi  which  excludes  terrain  conditions  will  be 
introduced  in  section  6.3,  and  then  a  new  terrain  quality  factor,  Qter-,  which  represents 
terrain  effects  will  be  introduced  in  section  6.4. 

6.3  LPI  Quality  Factors 

Using  the  expressions  for  the  communication  and  interception  ranges  given  previ¬ 
ously,  one  can  derive  an  expression  for  the  LPI  quality  factor,  QiPi,  as  follows: 


n2  _  PtGrcGcT ,  A  >2  1 

^  LcNsc  47r  Sc/Nsc 

(6.8) 

2  PiGpiGiT  ,  A  >2  1 

^  "  LrNsi  Si  in  si 

(6.9) 

(^C\2  _  ,GctGtc ^Nsi Si/Nsi  .,Li . 

Ri  GitG’I'i  Nsc  Sc/Nsc  Lc 

(6.10) 

The  LPI  quality  factor  is  defined  in  decibels  as 

Qr.Pi  =  =  20log&.  (6.11) 

IXJ  itj 

Also  the  LPI  quality  factor  can  be  written  as  the  sum  of  individual  subsystem  level  quality 
factors. 


20/og'(-— )  =  Qant  +  Qis  +  Qmod  +  Qatm 

Ki 

(6.12) 

where 

Qant  = 

in;  ,GctGtC\ 

Lr/T’CjTT’/ 

antenna  quality  factor 

(6.13) 

Qis  — 

interference  suppression  quality  factor 

(6.14) 
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modulation  quality  factor 


(6.15) 


^  107  r  ^ 

Qmod  — 

Qatm  =  lO^oS'(^)  atmospheric  quality  factor.  (6.16) 

Lc 

6.3.1  Antennas  Quality  Factor  (Qant)  ■  The  Qant  was  defined  as 

Qant  =  10/op(^F^)  (6.17) 

This  quality  factor  highlights  the  effect  of  the  antenna  gains  on  Q^pj.  If  the  Tx  uses  a 
directional  antenna  with  high  gain  in  the  desired  direction  and  small  side  lobes  or  nulls  in 
interceptor  direction,  then  a  large  Gtc  and  a  small  Gti  will  increase  Qant  and  reduce 
the  probability  of  interception.  This  forces  the  interception  Rx  is  move  closer  to  the  Tx  in 
order  to  maintain  the  original  probability  of  interception. 


6.3.2  Interference  Suppression  Quality  Factor  (Qis)-  The  Qjs  was  defined 
as 

Qjs  =  10/o5(^)  (6.18) 

This  quality  factor  is  basically  a  design  comparison  of  how  well  the  intended  Rx  suppresses 
interference  compared  to  the  interception  Rx.  If  iV5o,  the  total  interference  at  the  intended 
Rx  input,  is  lower  than  Nsj,  the  total  interference  at  the  interception  Rx  input,  then  the 
probability  of  interception  will  be  reduced  or  the  interception  Rx  is  forced  to  move  closer 
to  maintain  the  same  probability  of  interception. 


6.3.3  Modulation  Quality  Factor  {Qmod)  ■  The  Qmod  was  defined  as 

Qmod  =  l^log{  )  (6.19) 

This  is  the  most  important  factor  since  it  is  the  essence  of  system  performance.  In  the 
QmoD)  Sc  /  Nsc  is  the  SNR  required  for  a  specific  probability  of  bit  error.  Pi,,  in  the 


6-6 


intended  Rx,  whereas  5j  /  Nsi  is  the  SNR  required  for  the  same  Pb  in  the  interception 
Rx.  Usually,  the  interception  Rx  requires  a  lower  SNR  for  a  given  Pb-  However,  we  can 
increased  the  SNR  required  by  the  interception  Rx  and  decreased  the  SNR  required  by 
the  intended  Rx  by  using  special  techniques  such  as  transmitting  waveforms  with  large 
time-bandwidth  products  and  noise-like  spectra. 


6.3.4  Atmospheric  Quality  Factor  (Qatm)-  The  atmospheric  quality 

factor,  Qatmi  accounts  for  the  relative  atmospheric  effects  in  the  communication  and  in¬ 
tercept  links.  It  is  defined  as 


Qatm  —  —  Cc-Rc  (6.20) 

where  and  ^<7  are  loss  factors,  expressed  in  dB/km  or  dB/mile.  Sometimes,  the  path 
losses  in  the  intended  links  and  the  interception  links  can  be  assumed  to  be  equal  (ex¬ 
cluding  the  free  space  loss)  if  the  interception  and  the  intended  receiver  are  close  together 
and  operate  in  essentially  the  same  atmospheric  conditions.  In  these  cases,  Qatm  can 
be  neglected.  In  those  cases  where  Qatm  can  not  be  neglected,  the  reader  is  referred  to 
Ghordlo’s  work  [Ghordlo,  1996]  for  a  detailed  treatment  of  Qatm- 


6.3.5  Sample  LPI  Application  and  Limitations  of  Qlpi-  Solving  Eq.6.12 
for  Rj  yields 


20log{Ri)  =  2Qlog{Rc)  —  [Qant  +  Qis  +  Qmod  +  Qatm]  (6.21) 


or 


(Qaivt+<3;s+Qmod+Qatm]/20 


(6.22) 


_  yjQ(Q/lWT+0/S  +  QMOD+QATM]/20 

_ Rc _ 

2^Q(0.05)[QANT  +  Q/S+QMOD  +  OArM]  ’ 


(6.23) 

(6.24) 
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If  we  ignore  terrain  effects,  Ri  can  be  estimated  by  Eq.6.24.  In  this  case  Aj  is  a  circle 
centered  on  the  Tx.  Hence,  the  signal  can  be  intercepted  when  the  interception  Rx  is  inside 
the  circle  Aj.  Fig. 6. 3  shows  the  possible  interception  area  for  a  specific  transmitter  location 
and  three  different  values  of  Qant-  To  exemplify  the  analysis,  the  following  conditions 
were 

1.  Rc  is  the  maximum  possible  communication  range  and  is  fixed  at  5km. 

2.  To  compare  the  results,  three  different  Qant  (Oj  positive,  negative)  are  used. 

3.  The  other  quality  factors  are  assumed  to  be  0  (dB). 

Fig.6.3  demonstrates  how  one  can  reduce  Rj  (A/)  while  maintaining  the  original  Rq  {Ac) 
by  taking  advantage  of  the  relative  antenna  gains  in  the  pertinent  directions. 


Figure 

Qant 

Ri 

Note 

Fig.6.3.b 

0  (dB) 

5(km) 

Assume  Qis=Qatm=Qmod  =  O(dB) 

Fig.6.3.c 

lO(dB) 

1.6  (km) 

Fig.6.3.d 

-5(dB) 

8.9(km) 

Table  6.1  The  Qlpi  application  results. 


However,  this  form  of  Qipi  which  ignores  terrain  effects  is  not  sufficient  to  accurately 
characterize  the  ground  communication  because  terrain  effects  can  be  significant.  For  the 
current  Qlpi  definition,  changing  the  position  of  the  intended  Rx  along  an  arc  of  radius 
Rc  does  not  alter  Rj  or  Aj  .  This  is  demonstrated  in  Fig.6.4.  Note  that  the  results  are 
identical  to  those  shown  in  Fig.6.3. 

However,  if  terrain  effects  were  considered,  the  results  could  be  have  been  much 
different.  As  explained  in  chapter  5,  additional  path-loss  {L,hadow)  will  be  caused  by  the 
terrain  and  L shadow  is  very  sensitive  to  the  position  of  the  Tx  and  the  Rx.  Hence,  to  improve 
the  analysis  of  LPI  communication  over  rough  terrain,  the  terrain  effects  (Lghadow)  should 
be  considered.  In  the  next  section,  a  new  formulation  of  Qlpi  which  considers  terrain 
effects  will  be  introduced. 
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(a) 


Improbable 
Interception  Area 


Probable 
Interception  Area 


Figure  6.3  The  position  of  the  Tx  and  the  intended  Rx,  and  the  possible  interception 
area  (see  the  Table  6.1). 
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Distance  (inl 


The  Position  of  the  Tx  and  the  Intended  Rx 


Th«  tnt»nd»d R)r /a 
-  Tha  effffaranf  poaJ//cn 
•  Sama  Oafanoa  /SIrmJ 


y,  Intended  ' 


(a) 


Improbable 
Interception  Area 


Probable 
Interception  Area 


(b)Q>ijvT— O(dB) 


lO(dB) 


DltI«nM(m) 


(d)Q>lJVT — 5(dB) 


Figure  6.4 


The  different 


position  of  the  intended  Rx  and  the  possible  interception  area. 
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6-4  The  Terrain  Effect  and  Qlpi 

In  the  chapter  5,  UTCMOD  for  estimating  L,hadow  as  a  function  of  terrain  parameters 
{Mrei,  Vrei,  frequency,  etc.)  was  introduced.  In  the  past,  terrain  effects  were  ignored  in 
LPI  analysis  and  Lghadow  was  assumed  to  be  0  dB  in  all  directions.  As  shown  earlier,  the 
LPI  quality  factor,  Qlpi,  is  then  computed  as  follow: 

R 

Qlpi  =  20Zo5(^-)  =  Qant  +  Qis  +  Qmod  +  Qatm  (6.25) 

Ki 

However,  in  the  real  wireless  communication  scenarios,  L shadow  is  caused  by  signal  diffracted 
from  the  terrain  and  is  an  important  parameter  for  determining  communication  range. 
Thus  Qlpi  must  be  revised  to  include  a  terrain  quality  factor,  Qter,  which  accounts  for 
the  effects  of  terrain  on  signal  propagation  loss. 


LshadowidB)  =  M(dB)  +  V(dB)  +  R(dB) 


Link 

Condition 

M(dB) 

V(dB) 

R(dB) 

I 

(IVVerl)”-"®  •  (9.14  XlO-')  • 

-19  -f-lS.Zogiod 

11 /111 

(II^Terl)”-''  •  (5.33  XlO-®)  •  f 

-35  “I"  13*/o^iQci 

IV 

0 

0 

_ 

0 

Table  6.2  UTCMOD  diffraction  path-loss  parameters. 


6.4-1  LPI  Quality  Factor  {Qlpi)-  Including  terrain  effects,  the  received  power 
at  the  intended  Rx  {Sq)  is  estimated  as  follows: 


5c 


Pt^TcGcT^^ 

{^'nRcYLc[Lshadow]TC 


(6.26) 
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where 


•  Lc  is  the  atmospheric  loss  factor  {See  Fig. 6.1). 

•  [LahadowWc  is  the  diffraction  path-loss  between  the  Tx  and  the  intended  Rx. 
Also  the  received  signal  power  to  noise  PSD  is 

Sc  _  PiGtcGct  /  ^2 

NsC  LcNsc[Lshadow\TC  47ri?c 

Prom  Eq.6.27,  solving  for  the  maximum  communication  range  (i?c)  yields 


(6.27) 


I  PiGtcGct 

A2  1  1 

y  LcNsc 

(47r)2  Sc/Nsc  [Lahadowlrc 

(6.28) 


Ri  can  be  estimated  with  similar  steps.  The  received  power  at  the  interception  Rx  is 

PtGriGjT^^ 


Si 


{AtxRiYLi  [Lahadow]  TI 


(6.29) 


where 


•  Lj  is  the  atmospheric  loss  factor  between  the  Tx  and  the  intercept  Rx  (^ee  Fig. 6.1). 

•  [L>ahadow\Ti  is  the  diffraction  path-loss  between  the  Tx  and  the  interception  Rx. 

The  received  signal  power  to  noise  PSD  can  be  expressed  in  the  terms  of  the  link  parame¬ 


ters: 


PiGtiGit  /  a  2 


Nsi  LiNsi[Lahadow\TI  47ri?/' 

From  Eq.6.30,  solving  for  the  maximum  interception  range  {Ri)  yields 


(6.30) 


Ri 


!  PiGtiGit  A^ 


LiNsi  {^T^fSjlNsi 


[L 


$hadow\TI 


Dividing  Eq.6.28  by  Eq.6.31  and  squaring  the  results  yields 


,Pc  ',2  _  ,GctGtC  Y  Nsi  w  Si/Nsi  ..  Li  [Lahadow]TI  >. 

Ri  GitGti  Nsc  Sq/Nsc  Lc  [Lahadow\TC 


(6.31) 


(6.32) 
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Finally,  the  LPI  quality  factor  which  includes  the  terrain  effects  is  defined  as  follows: 


Qlpi  —  Qant  +  Qis  +  Qmod  +  Qatm  +  Qter  (6.33) 

where 

)  (6.34) 

Y^shadowlTC 

6.4-2  Terrain  Quality  Factor  (Qter)-  In  rough  terrain,  the  terrain  statistics 
are  sensitive  to  the  position  of  the  Tx  and  the  Rx.  Hence,  the  terrain  statistics  between 
the  Tx  and  the  intended  Rx  are  different  from  the  terrain  statistics  between  the  Tx  and 
the  interception  Rx.  Also,  the  specific  link  condition  between  the  Tx  and  the  Rx  influences 


Figure  6.5  The  terrain  parameters  between  the  Tx  and  the  Rx 


■^shadow 


(see  Table  6.2).  Thus,  Qter  can  be  estimated  by  the  following  steps: 


1.  Determine  the  positions  of  the  Tx,  the  intended  Rx,  and  the  interception  Rx. 

2.  Estimate 


•  Determine  the  terrain  statistics  between  the  Tx  and  the  intended  Rx  (Fig.6.5). 


•  Check  the  link  condition  between  the  Tx  and  the  intended  Rx. 
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•  Estimate  the  [Lahadow\TC  with  Table  6.2. 

3.  Estimate 

•  Determine  the  terrain  statistics  between  the  Tx  and  the  interception  Rx  (Fig.6.5). 

•  Check  the  link  condition  between  the  Tx  and  the  interception  Rx. 

•  Estimate  the  [Lshadow\Ti  with  Table  6.2. 

4.  Estimate  the  Qter- 

Qter  =  (6.35) 

\:^8hadow\TC 

The  terrain  quality  factor  (Qter)  is  sensitive  to  the  positions  of  the  Tx,  the  intended  Rx, 
and  the  interception  Rx.  The  Fig.6.6.a  shows  the  position  of  the  Tx  and  the  intended  Rx. 
Fig.6.6.b  shows  the  height  variation  of  the  terrain.  For  fixed  Tx  and  intended  Rx  locations, 
Qter  was  computed  as  a  function  of  the  position  of  the  interception  Rx.  Fig. 6. 7  shows 
Qter  for  each  specific  point  assuming  the  interception  Rx  is  on  that  point.  Since 

‘20log{-^)  =  Qant  +  Qis  +  Qmod  +  Qatm  +  Qter-  (6.36) 

Then,  one  can  calculate  Rj  for  a  fixed  Tx  position  and  each  fixed  intended  Rx  position. 
Then  Aj  can  be  computed  from  iZ/.  However,  Aj  is  no  longer  a  valid  interception  metric 
because  the  Qter  (or  [L^hadowWc  and  [LghadowWi)  are  very  sensitive  to  the  position  of  the 
Tx,  the  intended  Rx,  and  the  interception  Rx.  Therefore,  the  probability  of  interception 
is  estimated  not  by  the  range  (or  the  area)  but  by  the  relative  positions  of  the  Tx,  the 
intended  Rx,  and  the  interception  Rx.  In  the  next  section,  a  new  interception  metric  will 
be  introduced. 

6.5  Tactical  Quality  Factor  (Qtac) 

LPI  communication  theory  has  been  developed  to  characterize  the  military  commu¬ 
nication  objective  of  maximizing  Rc  with  respect  to  Rj.  Using  the  existing  Qlpi  which 
ignores  terrain  conditions  (roughness),  Rj  can  be  computed  with  Eq.6.21  through  Eq.6.24. 
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Specifically 


(6.37) 


^  '[Qi^-^^)lQANT  +  QlS+QMOD~\-QATM] 

The  interception  Rx  has  a  high  probability  of  detecting  the  signal  when  the  range  between 
the  Tx  and  the  interception  Rx  is  less  than  Rj. 

However,  in  the  real  communication,  this  is  not  an  accurate  characterization  because 
L shadow  and  hence  Qter  are  very  sensitive  to  the  relative  transmitter  and  receiver  location. 
Hence,  a  new  measure  which  estimates  the  probability  of  the  interception  is  desired.  In 
this  thesis,  the  tactical  quality  factor  (Qtac)  is  suggested  for  this  new  measure. 


6.5.1  Basic  Theory  of  Qtac-  Eq.6.27  represents  the  required  signal-to- noise 
ratio  for  a  given  probability  of  detection  at  the  intended  Rx.  Prom  Eq.6.27,  the  required 
transmit  power,  Pt,  can  be  estimated. 


{Pt)sc  = 


Sc  Lc[Lshadow]TcNsC  ^2 
NsC  Gtc^cT  ^ 


(6.38) 


where  {Pt)sc  is  the  required  transmit  power  from  the  Tx  to  the  intended  Rx. 
Similarly,  the  Pt  required  by  the  interception  Rx  is 


{Pi)si  = 


Sj  Li[Lghaaow]TjNsi 
Ns  I  GtiGjt 


4^  2 

^  A  ’ 


(6.39) 


where  {Pt)si  is  the  required  transmit  power  from  the  Tx  to  the  interception  Rx. 
Then  the  ratio  {Pt)si  to  {Pt)sc  can  be  written  as  follows: 


{Pi)si 

{Pt)sc 


S,  Lj\L  /  477 ilf  \2 

Nsi  GtiG  iT  ^  A  ' 

Sc  Lc{Lahad  ou>  Irc^sc  (A-nRc  \2 

Nsc  GtcGct  ^  a  ' 

,GctGtc Nsi  w  Sj/Nsi  ^ .  Lj  . ,  {Lshadow]Ti  w  Ri  ',2 
GjtGti  Nsc  Sc/Nsc  Lc  [L  shadow]TC  Pc 


(6.40) 

(6.41) 


In  Eq.6.41,  the  ratio  {Pt)si  to  {Pt)sc  is  a  function  of  the  quality  factors  and  the  distances. 
Recall  that  L shadow  is  also  a  function  of  distance,  i?/  is  the  distance  between  the  Tx  and  the 
interception  Rx  whereas  Rc  is  the  distance  between  the  Tx  and  the  intended  Rx.  Since, 
is  a  result  of  the  spherical  wave  spreading  out  or  what  is  commonly  called  free  space 
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loss,  I  will  call  this  the  free  space  distance  quality  factor  (Qdis)-  Qdis  in  dB  given  by: 

Qdis  =  20Zop(-^).  (6.42) 

(6.43) 

where 

Qsys  —  Qant  +  Qis  +  Qmod  (6.44) 

and 

Qenv  =  Qatm  +  Qter  +  Qdis  (6.45) 

Note  that  Qsys  is  exactly  the  old  Qlpi  and  only  contains  system  design  parameters. 
Qenv  I  on  the  other  hand,  contains  all  of  the  parameters  that  depend  on  the  particular 
operating  environment  of  the  system.  Qtac  represents  the  analytical  relationship  between 
{Pi)sc  and  (Pt)si-  If  tho  Qtac  is  10  dB,  then  the  interception  Rx  cannot  intercept  the 
signal  on  that  point  unless  the  Tx  increases  the  transmitted  power  lOdB.  On  the  other 
hand,  if  Qtac  is  -10  dB,  then  the  interception  Rx  can  intercept  the  signal  at  that  point  until 
Tx  decreases  the  transmitted  power  lOdB.  However,  after  decreasing  the  transmitted  power 
lOdB,  the  intended  Rx  also  cannot  detect  the  signal.  Thus,  friendly  communications  must 
change  their  tactics  (e.g.  move  closer  together)  to  obtain  a  positive  Qtac  and  deny  the 
enemy  an  opportunity  to  intercept  the  friendly  signal.  Hence,  Qtac  is  an  new  approach 
for  analyzing  tactical  communication  (COMM)  and  electronic  warfare  (EW)  scenarios 
that  accounts  for  the  effects  of  terrain  and  possibly  allows  one  to  exploit  terrain  to  create 
a  COMM/EW  advantage.  In  the  next  section,  examples  of  Qtac  applications  will  be 
introduced. 

6.5.2  Sample  Qtac  Applications  for  Analyzing  the  LPI  Performance  of 
a  Tactical  Communication  System.  Qtac  represents  the  analytical  relationship 
between  the  required  {Pt)sc  and  iPt)si  based  on  terrain  effects  and  all  of  the  other  qual- 


Finally,  the  tactical  quality  factor  is  defined  as 


Qtac  —  'i-0log( 


{Pt)si 

{Pt)sc 


)  —  Qsys  +  Qenv 
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ity  factors.  Since  Qtac  is  sensitive  to  the  positions  of  the  Tx  and  the  Rx,  it  has  wide 
applicability  for  analyzing  tactical  communication  (COMM)  and  electronic  warfare  (EW) 
scenarios  for  the  wireless  communication.  In  this  thesis,  an  example  of  how  to  estimate 
potential  interception  positions  (areas)  and  how  to  optimize  the  Tx  position  will  demon¬ 
strate  the  power  of  this  approach  for  analyzing  LPI  communication  systems  and  tactical 
implementations. 


6. 5. 2.1  Possible  Interception  Regions.  In  the  section  6.3.5,  the  po¬ 

tential  interception  area  was  estimated  with  the  previous  Qlpi  definition  which  ignores 
terrain  effects.  In  this  section,  the  potential  interception  area  will  be  estimated  with  Qtac 
and  compared  to  the  results  in  the  section  6.3.5  (Fig.6.3). 

For  fixed  Tx  and  intended  Rx  positions  (Fig.6.8.a),  Qtac  for  each  point  can  be  read¬ 
ily  calculated  and  displayed  (Fig.6.8.b).  At  points  where  Qtac  is  greater  than  20  dB,  the 
interception  Rx  cannot  detect  the  signal  unless  the  Tx  power  increases  by  at  least  20  dB. 
Hence,  the  interception  Rx  must  do  something  (e.g.  increase  antenna  gain)  to  increase  the 
performance  of  the  intercept  link  or  move  to  more  favorable  location  (i.e.  Qtac  <  0)  to  de¬ 
tect  and  intercept  the  signal.  However,  at  the  point  where  Qtac  is  -30  dB,  the  interception 
Rx  detects  the  signal  until  the  Tx  power  decreases  by  at  least  30  dB.  Of  course,  this  loss  in 
transmit  power,  Pt,  must  be  compensated  for  by  a  comparable  gain  in  some  other  quality 
factor  to  ensure  successful  communication  on  the  intended  link.  Therefore,  at  any  points 
where  Qtac  is  less  than  0  dB  the  interception  Rx  detects  the  signal  and  these  points  com¬ 
prise  the  potential  interception  region.  Fig.6.8.c  shows  the  most  likely  interception  regions 
predicted  by  Qtac-  Fig.6.8.d  shows  the  estimated  potential  interception  area  predicted  by 
Qlpi-  For  simplicity  in  this  example,  only  Qter  and  Qdis  are  considered  while  the  other 
quality  factors  are  assumed  to  be  0  dB.  Clearly,  the  terrain  effects  are  apparent  in  Fig.6.8.c 
whereas  they  are  absent  in  Fig.6.8.d.  This  strikingly  demonstrates  the  superiority  of  using 
the  tactical  quality  factor,  Qtac-,  as  an  LPI  performance  assessment  tool  compared  to  the 
previously  accepted  LPI  quality  factor,  Qlpi,  approach. 
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(c) 


Figure  6.8 


DI»ltnM(tn) 


Most  likely  interception  region 
predicted  by  Qtac 


(d).  Potential  interception  area 
predicted  by  Qlpi 


Improbable 
Interception  Area 


Probable 
Interception  Area 


A  comparison  of  Qtac  and  Qipj  predictions  of  region  vulnerable  to 
interception. 


6-19 


6. 5. 2.2  Using  Qtac  to  Locate  Optimum  LPI  Transmitter  Positions. 

In  a  tactical  communication  scenario,  finding  the  Tx  position  that  affords  an  advantage 
to  the  intended  Rx  and  a  disadvantage  to  the  interception  Rx  would  be  very  useful  to 
strategists,  tacticians,  and  communication  system  planners.  If  we  know  the  interception  Rx 
position,  then  advantageous  transmitter  locations  can  be  determined  using  the  Qtac  LPI 
metric.  Fig.6.9  shows  Qtac  and  the  suggested  Tx  positions  for  the  anticipated  interception 
and  intended  Rx  positions.  For  simplicity  in  this  example,  only  Qter  and  Qdis  were 
considered  and  while  other  quality  factors  were  assumed  to  be  0  dB.  Again,  the  power  of 
the  Qtac  metric  is  readily  apparent  as  the  results  depicted  in  Fig.6.9.c  could  be  a  valuable 
addition  to  a  commander’s  assessment  of  the  battle  field. 


The  position  of  the  Intended  Rx  and  the  Interception  Rx 
1km  a''' 


f 
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of  the  j 
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(a). The  Rx  positions 


(b).QxAC  for  every  possible  Tx  location 


Taotloal  Quality  Factor 


strongly  Suggested  Point  (Qtac  >  20  dB) 


Suggested  Point  (0  dB  <  Qtac  <  20dB ) 


Dissuaded  Point  (-20  dB  <  Qtac  <  0  dB) 


Strongly  Dissuaded  Point  (Qtac  <  -20  dB) 
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(c)  Final  Tx  position  recommendation. 


Figure  6.9  Tx  location  recommendations  using  Qtac- 
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6.6  Summary 

In  military  wireless  communication  scenarios,  terrain  can  be  used  to  create  a  com¬ 
munication  advantage  for  friendly  forces.  By  using  a  newly  developed  stochastic  model 
(UTCMOD)  which  estimates  the  terrain  effects  (Lshadow)  on  radiowave  propagation,  many 
tactical  communication  (COMM)  and  electronic  warfare  (EW)  scenarios  can  be  readily 
analyzed  and  optimized  for  LPI  performance.  Specifically,  by  maximizing  Qteri  we  can 
use  terrain  to  achieve  a  tactical  communication  advantage. 

Previous  Qipi  definitions  are  overly  simplistic,  too  sensitive  to  range,  and  ignore 
terrain  effects.  Qtac  on  the  other  hand,  overcomes  the  Qlpi  shortcomings.  Two  Qtac 
applications  were  demonstrated  to  exhibit  the  tactical  utility  of  this  metric.  Using  Qter 
and  Qtac  to  assess  LPI  performance  allows  one  to  exploit  terrain  to  achieve  a  tactical  com¬ 
munication  advantage  by  creating  sub-optimum  conditions  for  the  enemy  and  ensuring  the 
best  possible  conditions  for  friendly  forces. 
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VII.  Conclusions  and  Recommendations 


This  research  originated  with  the  intent  to  analyze  terrain  effects  on  signal  propaga¬ 
tion  and  to  improve  the  analysis  of  tactical  communication  scenarios  over  rough  terrain. 
Thus  the  objectives  of  this  research  were  the  following: 

1.  Apply  a  reasonable  propagation  model  for  analyzing  the  mean  path-loss  of  signal 
propagation  over  rough  mountainous  terrain  that  is  valid  over  various  frequency 
ranges. 

2.  Provide  a  valuable  tool  and  insight  for  analyzing  LPI  characteristics  of  ground  com¬ 
munication  systems. 

These  objectives  were  achieved  and  the  results  of  this  research  exceeded  expectations.  This 
chapter  summarizes  the  results  and  conclusions  and  recommends  future  research  motivated 
by  what  was  discovered  during  this  effort. 


7.1  Summary 

Universal  Terrain  Channel  Model  (UTCMOD),  a  broad  banded  propagation  model 
for  estimating  the  mean  path-loss  of  signal  propagation  over  rough  mountainous  terrain, 
was  developed  by  simulating  terrain  effects  on  signal  propagation. 

Using  UTCMOD,  many  tactical  communication  (TAC  COMM)  and  electronic  war¬ 
fare  (EW)  scenarios  can  be  analyzed.  Furthermore,  two  new  quality  factors,  Qter  and 
Qtacj  resolved  numerous  shortcomings  of  a  previously  defined  LPI  quality  factor  or  Qlpi 
This  new  method  of  LPI  analysis  allows  one  to  exploit  terrain  to  create  a  communication 
advantage  for  friendly  forces. 


7.1.1  UTCMOD  Summary.  The  prediction  of  path-loss  is  a  very  important 
step  in  planning  an  effective  ground  wireless  radio  system.  Hence,  several  propagation 
models  which  are  based  on  experimental  results  have  been  developed  for  point-to-point 
mobile  communication.  Of  these  models,  the  Hata  model,  a  widely  used  model,  was  in- 
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troduced  in  Chapter  11.  However,  the  Hata  model  is  not  sufficient  for  predicting  path-loss 
in  rough  terrain.  The  model  is  most  relevant  for  terrain  conditions  characteristic  of  large 
cities  or  medium-small  cities.  Thus  it  is  inadequate  for  many  different  kinds  of  terrain  con¬ 
ditions  that  may  be  encountered.  Hence,  a  channel  model  that  estimates  mean  path-loss 
over  many  different  kinds  of  terrain  conditions  is  desired. 

To  derive  a  reasonable  propagation  model,  Etotai  fhc  receiving  point  based  on 
McNamara’s  signal  diffraction  theory  was  discussed  in  Chapter  III.  After  comparing  com¬ 
putations  using  Balanis’  wedge  equation  with  computations  using  McNamara’s  knife-edge 
equation,  it  was  shown  that  the  McNamara  knife-edge  equation  is  reasonable  for  a  variety 
of  wedge-like  obstacles.  Using  McNamara  diffraction  calculations.  Real  Terrain  Diffraction 
Model  (RTDMOD)  was  developed.  RTDMOD  can  predict  path-loss  from  real  terrain  data 
and  is  not  limited  by  the  number  of  obstruction.  However,  RTDMOD  is  computationally 
intensive.  Hence,  a  computationally  efficient  stochastic  propagation  channel  model  was 
created  that  is  easy  to  apply  to  any  terrain  conditions. 

UTCMOD,  developed  using  RTDMOD,  is  a  reasonable  propagation  model  for  esti¬ 
mating  the  mean  path-loss  between  the  Tx  and  the  Rx  and  it’s  validity  was  tentatively 
verified  by  comparing  the  results  with  the  Hata  model.  To  estimate  the  terrain  effects  on 
signal  propagation,  UTCMOD  has  the  following  desirable  qualities: 

1.  UTCMOD  is  applicable  to  many  different  kinds  of  terrain  conditions.  The  model  es¬ 
timates  the  additional  path- loss  from  the  mean  and  the  variance  of  the  terrain  height, 
the  propagation  distance,  and  the  carrier  frequency  (Table  5.8).  Since  the  terrain 
statistics  (mean  and  variance)  determine  the  path-loss,  UTCMOD  is  applicable  to 
many  different  kinds  of  terrain  conditions  whereas  the  Hata  model  is  most  relevant 
for  terrain  conditions  characteristic  of  large  cities  or  medium-small  cities. 

2.  UTCMOD  considers  the  communication  link  condition.  In  addition  to  the  fundamen¬ 
tal  definitions  of  communication  modes  (Table  5.1),  the  communication  link  condi¬ 
tions  are  important  for  estimating  the  path-loss.  New  expanded  communication  link 
definitions  described  in  chapter  5  (Table  5.2)  provide  the  correct  perspective  for 
analyzing  the  propagation  effects. 
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3.  UTCMOD  is  valid  over  various  frequency  ranges  and  estimates  the  mean  path-loss 
of  spread  spectrum  (SS)  signal  propagation.  SS  signals  have  a  frequency  range  (/; 
<  fp  <  fh)-  If  we  assume  the  probability  of  occurance  for  each  frequency  between 
fi  and  fh  is  the  same,  then  the  path-loss  can  be  approximated  by  UTCMOD  (Table 
5.11). 

7.1.2  Tactical  Communication  (TAG  COMM)  and  Electronic  Warfare 
(EW)  Scenarios.  In  military  tactical  communication  (TAG  COMM)  and  electronic 
warfare  (EW)  scenarios,  Low  Probability  of  Interception  (LPI)  is  an  important  issue.  Thus, 
LPI  communication  theory  has  been  developed  throughout  the  past  decade.  However,  the 
existing  LPI  theory  is  not  sufficient  for  real  wireless  communication  because 

1.  The  previous  analysis  methods  using  Qlpi  ignore  the  terrain  effects  on  signal  prop¬ 
agation.  Thus,  if  the  position  of  the  intended  Rx  is  changed  but  Rq  and  the  other 
quality  factors  are  held  constant,  then  the  predicted  interception  area  will  be  the 
same.  However,  in  real  wireless  communication  that  is  not  true  because  variable 
terrain  effects  will  change  the  potential  interception  regions. 

2.  Qipr,  as  previously  defined,  is  incompatible  with  other  quality  factors  which  are  a 
function  of  distance.  The  current  Qlpi  assesses  LPI  performance  using  the  (^) 
metric.  Thus,  if  any  quality  factor  is  a  function  of  distance,  the  existing  Qlpi  metric 
does  not  work.  As  a  example,  if  Qatm  =  ^log{Rj)  -  3log{Rc),  then 

2Qlog{Rc)  —  20log{Ri)  =  Qant  +  Qis  +  Qmod  +  ‘llog{Rj)  —  3log{Rc)  (7.1) 
23log{Rc)  —  2^log{Ri)  =  Q ant Qis  Qmod-  (7-2) 

Now,  it  is  very  difficult  to  estimate  the  LPI  performance  from  the  Eq.7.2.  Further¬ 
more,  if  a  new  quality  factor  is  introduced  that  is  explicitly  a  function  of  distance, 
then  it  is  very  difficult  to  estimate  the  LPI  performance  by  adding  the  new  quality 
factor  to  the  existing  Qlpi- 

Fortunately,  these  problems  can  be  resolved  by  using  Qter  and  Qtac-  Qter  represents 
the  path- loss  relationship  between  [LghadowlrCi  the  additional  diffraction  path- losses  be- 
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tween  the  Tx  and  the  intended  Rx,  and  [Lahadowlri,  the  additional  diffraction  path-loss 
between  the  Tx  and  the  interception  Rx.  Qtac  represents  the  analytical  relationship  be¬ 
tween  {Pt)sci  the  required  transmitted  power  for  the  communication  link,  and  {Pt)sij  the 
required  transmitted  power  for  the  interception  link.  Qtac  was  based  on  the  existing  LPI 
performance  metric,  Qlpi,  and  Qter-  Hence,  the  first  problem  associated  with  Qlpi  can 
be  solved  with  Qter  and  the  second  problem  with  Qlpi  can  be  solved  with  Qtac-  Fi¬ 
nally,  Qter  and  Qtac  can  be  used  to  analyze  tactical  communication  (TAG  COMM)  and 
electronic  warfare  (EW)  scenarios  in  order  to  exploit  terrain  to  create  a  communication 
advantage  for  friendly  forces. 

Again  it  is  worth  emphasizing  that  although  this  work  caste  in  a  military  perspective, 
all  of  the  models,  quality  factors,  and  concepts  are  directly  applicable  in  currently  expand¬ 
ing  commercial  applications  involving  signal  propagation  in  mountainous  rural  areas. 


7.2  Recommendations 

This  thesis  provided  a  stochastic  propagation  channel  model  (UTCMOD)  and  high¬ 
lighted  it’s  usefulness  in  tactical  situations.  Throughout  this  thesis,  only  the  two  dimen¬ 
sional  (2-D)  signal  diffraction  (forward  diffraction)  was  considered  for  estimating  the  effects 
of  terrain.  In  2-D  diffraction  theory,  only  the  elevation  angles  between  the  Tx,  obstacle, 
and  the  Rx  are  considered  when  calculating  the  total  field  at  the  receiving  point  (Efotai)- 
However,  the  signals  are  omni-directionally  diffracted  at  the  diffraction  point  (obstacle) 
and  the  Rx  can  detect  the  sideward  diffracted  signals  that  was  diffracted  by  an  obstacle 
at  an  azimuth  angle  that  is  different  from  the  azimuth  angle  of  the  Tx.  If  we  consider 
very  rough  ground  conditions,  sometimes  the  path-loss  of  the  sideward  diffracted  signals 
can  be  lower  than  the  path-loss  of  the  forward  diffracted  signal.  Hence,  expanding  the 
analysis  to  include  three  dimensional  (3-D)  diffraction  theory  is  recommended  to  improve 
the  performance  assessment  of  ground  wireless  communications. 

In  3-D  diffraction  theory,  both  of  the  elevation  and  the  azimuth  angle  are  considered 
to  calculate  the  total  field  at  the  receiving  point  {Etotai)-  Also  the  3-D  diffraction  channel 
model  should  focus  on  finding  the  elevation  and  the  azimuth  angle  of  the  diffraction  point 
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which  causes  the  lowest  propagation  path-loss.  If  we  can  determine  these  critical  angles, 
then  the  performance  of  the  radio  system  can  be  greatly  improved  by  adjusting  the  gain 
of  the  antenna  in  that  direction. 
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Appendix  A.  Abbreviations  and  Symbols 
A.l  Abbreviations 

CDMA:  Code  Division  Multiple  Access 

DG  Method:  Daygout  Method 

DSSS:  Direct  Sequence  Spread  Spectrum 

EC:  Electronic  Combat 

EP  Method:  Epstein-Peterson  Method 

EW:  Electronic  Warfare 

FHSS:  Frequency  Hopping  Spread  Spectrum 

FSK:  Frequency  Shift  Keying 

GTD:  Geometric  Theory  of  Diffraction 

LOS:  Line  of  Sight 

LPI:  Low  Probability  of  Interception 

MBPLG:  Mean  based  path-loss  generator 

MVBPLG:  Mean-variance  based  path-loss  generator 

MVWBPLG:  Mean-variance-wavelength  based  path-loss  generator 

MVWDBPLG:  Mean-variance-wavelength-distance  based  path-loss  genera¬ 
tor 

PEC:  Perfect  Electrically  Conducting 
PN:  Pseudo-Noise 

RTDMOD:  Real  Terrain  Diffraction  Model 

RTG:  Random  Terrain  Generator 

SS:  Spread  Spectrum 

TAC  COMM:  Tactical  Communication 
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THSS:  Time  Hopping  Spread  Spectrum 
UTCMOD:  Universal  Terrain  Channel  Model 

A.  2  Symbols 


a{hre,  fc)  '•  Correction  factor  in  Hata  model 
A:  Spatial  attenuation 
Ar'.  Effective  area  of  the  receiver  antenna 
/3:  Wave  number  (—  2  tt  /  X) 

C(x):  Fresnel  cosine  integral 
V:  Diffraction  coefficient. 

d:  Propagation  distance  between  transmitter  and  receiver 

A/:  Hop  frequency  spacing 

^dif fraction*  Diffraction  field 

E incidence-  Incidence  field 

Etotai-  Total  field  at  receiving  point 

{Etotai)A-  Total  field  at  the  receiving  point  in  actual  communication 
{Etotai) F-  Total  field  at  the  receiving  point  in  free  space 
Fi'.  First  Fresnel  zone  radius 
/<,:  Carrier  frequency 

fh'.  Highest  frequency  of  the  propagation  bandwidth 
/;:  Lowest  frequency  of  the  propagation  bandwidth 
Gr'.  Receiver  antenna  gain 
Gt'.  Transmitter  antenna  gain 
he'.  Path  clearance 
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hre-  Receiver  antenna  height 
hte-  Transmitter  antenna  height 
Htx  ■  Terrain  height  of  the  transmitter 
A:  Signal  wavelength(= 

L atmosphere-  Atmospheric  path- loss 
Lfree-  Free  space  path- loss 
Liii'.  Lit  diffraction  path- loss 
L shadow-  Shadow  diffraction  path-loss 

[L3hadow]TC-  Shadow  diffraction  path-loss  between  the  transmitter  and  the  in¬ 
tended  receiver 

[LshadowjTi-  Shadow  diffraction  path-loss  between  the  transmitter  and  the  in¬ 
terception  receiver 

LSS shadow-  Shadow  diffraction  path-loss  for  Spread  Spectrum 

Mter-  Mean  of  terrain  height 

Mtx-  Relative  mean  of  terrain  height 

4>'-  Observation  angle.  Angle  between  obstacle  and  the  receiver 

0’:  Incidence  angle.  Angle  between  the  transmitter  and  obstacle 

Pav-  Average  received  signal  power  density 

Pb  '■  Probability  of  bit  error 

Pr'.  Received  signal  power 

Pf.  Transmitted  signal  power 

Qant-  Antenna  quality  factor 

Qatm-  Atmospheric  quality  factor 

Qdis-  Free  space  distance  quality  factor 

Qenv-  Environmental  quality  factor 
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Qjs'.  Interference  suppression  quality  factor 

Qlpi'  LPI  quality  factor 

Qmod-  Modulation  quality  factor 

Qsys-  System  quality  factor 

Qtac-  Tactical  LPI  quality  factor 

Qter-  Terrain  LPI  quality  factor 

Re'.  Chip  rate 

Re'.  Communication  range 

Ri'.  Interception  range 

Rx:  Receiver 

S(x):  EYesnel  sine  integral 

Tx:  Transmitter 

v:  Diffraction  parameter 

Vp'.  Velocity  of  propagation 

Vter  ■  Variance  of  terrain  height 

Wd'.  Signal  bandwidth  before  spreading 

Wss'.  Signal  bandwidth  after  spreading 
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Appendix  B.  KRAUS  Diffraction  Method 

Huygens’  principle  of  physical  optics  can  be  used  to  explain  the  apparent  bending 
of  radio  waves  aronnd  obstacles, i.e.,  the  diflFraction  of  waves.  A  diffracted  ray  is  one  that 
follows  a  path  that  cannot  be  interpreted  as  either  reflection  or  refraction  [Kraus,  1992]. 
As  an  example,  consider  a  uniform  plane  wave  incident  on  a  conducting  half-plane. 


_ _ !L  I  " _ '[ _ "  |  " 

Conducting  P'®"®  *®''® 


half-plane 


pteliiii 


Illuminated 

side 


ib) 


^Physical 

optics 


Ka  (k  •  constant) 


Figure  B.l  Plane  wave  incident  from  above  onto  a  conducting  half-plain  with  resultant 
power-density  variation  below  the  plane  as  obtained  by  physical  optics. 


as  in  Fig.B.l.  We  want  to  calculate  the  electric  field  at  point  P  by  using  Huygens’  principle; 


E=  f  dE 

Jover  X  axis 


(B.l) 


where  dE  is  the  electric  field  at  P  due  to  a  point  source  at  a  distance  x  from  O,  as  in 


where  r  is  the  distance  from  obstacle  (m)  and  a  is  the  distance  into  shadow  region  (m).  If 
6  T,  follows  that 

6  =  ^  (BA) 


2r 


When  we  let  P  =  2/rA  and  u  =  kx,  Eq.B.3  becomes 


JP  /•OO 

E  =  l^du  (B.5) 

kr  Jka 

which  can  be  rewritten  as 

TP 

E  =  (  /  e-^"“  '^du  -  /  6-^’^“  l^du)  (B.6) 

kr  Jo  Jo 

The  Integrals  in  Eq.B.6  have  the  form  of  Fresnel  integrals  [A  6ramowit2:  and  Stegun,  1964]. 


Cix)  = 

-KV?  , 

/  cos-—-du 

Jo  2 

Fresnel  cosine  integral 

(B.7) 

S(x)  = 

r  .  7^2 

/  sm——du 

Jo  2 

Fresnel  sine  integral 

(B.8) 

Also,  from  the  Fig.B.2,  C(oo)  «  S(oo)  «  So  that  Eq.B.6  can  be  written 


Fresn^  Integrals 


Figure  B.2  Fresnel  Integrals. 


^  ( I  +  >i  -  +jS{ka)] ) 


(B.9) 
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where 


Appendix  C.  Balanis’  Wedge  Diffraction  Equation 


In  his  reference  [Balanis^  1989],  Balanis  represents  total  field  at  the  receiving  point 
in  the  following  form: 


Etotai  =  E^-V-A-e-^^^ 


(C.l) 


where 


Eo  =  Free  space  field  strength 
V  =  Diffraction  coefficient 
A  =  Spatial  attenuation 
p  =  2-k/X 

s  =  Propagation  distance,  m 


C.l  Diffraction  Coefficient  (V) 


observation 


Figure  C.l  Incidence  angle(^’)  and  Observation  angle((/>). 


C-1 


Fig.C.l  shows  the  incidence  angle((^’)  and  the  observation  angle(^).  The  diffraction 
coefBcient(r>)  can  be  calculated  with  the  incidence  angle(0’)  and  the  observation  ang\e{(j)). 

X>(p,  (f),  (f)',  n)  =  T>\p,  (f)  -  <f>',  n)  -  'D''{p,  (j)  +  </>',  n)  (C.2) 


where  angle  of  wedge  is  (2-n)7r  and  n  can  be  any  real  number.  Also 


V\p,(l>- (j)\n) 

(C.3) 

(C.4) 

V''{p,(t>  +  4>\n) 

■  2„vS3'“’‘'  \-nM{4'  +  4>)] 

(C.5) 

(C.6) 

where 

=  1  +  cos\{(j)  ±  (f))] 

(C.7) 

T[i5pg{ct>^4>')\ 

=  /  e-^^"du 

J  y/0pG{<t>±<t>') 

(C.8) 

Using  the  computer  program,  the  diffraction  coefficient('P)  can  be  calculated. 

However, 

the  diffraction  coefficient (X>)  is  a  function  of  cot[7v  ±  ±  (f>')].  Thus,  ii  (j)  ±  (f)'  =  mvr  for 

any  integer  m,  then  cot[K  ±{(f)±  is  infinite  and  the  computation  is  invalid. 


C.2  Spatial  Attenuation  (A) 

Cylindrical  wave  incidence  is  considered  for  calculating  spatial  attenuation(v4).  Thus 


A  = 


1 


(C.9) 


where  p  is  the  distance  between  the  diffraction  point  and  the  observation  point  and  /?o  is 
the  angle  as  shown  in  Fig.C.2. 
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Figure  C.2  Distance  and  angle  for  calculating  the  spatial  attenuation 


C-3 


Appendix  D.  Lit  Diffraction  {Lm)  for  Multiple  Subpath  Edges 
D.l  Basic  Theory  for  Computing  Lm 

Eq.3.2  represents  the  total  electric  field  at  the  receiving  point.  In  the  lit  region,  the 
total  field  is  composed  of  an  incidence  field  and  a  diffraction  field.  The  diffraction  field 
due  to  one  subpath  obstacle  can  be  estimated  by  subtracting  the  incidence  field  from  the 
total  field. 


^dif  fraction  — 


total 


-E, 


incidence 


(D.l) 


The  total  field  in  the  lit  region  due  to  multiple  subpath  obstacles  is 

N 

Efotal  —  ^  ]  Ej  Emcidence 
1=1 

where  N  is  the  number  of  lit  edges  (subpath  obstacles)  and  Ei  is  the  diffraction  field, 
Ediffraction,  due  to  the  i-th  lit  edge  (subpath  obstacle).  In  Fig.D.l,  the  total  field  at  the 
Rx  is  composed  of  an  incidence  field  and  three  different  diffraction  fields. 

However,  Eq.D.2  is  not  practical  for  calculating  Lm  over  rough  terrain  because  of  the 


Total  Field  at  the  Rx  for  Multiple  Subpath  Edges 


Tx  Edge  1  Edge  2  Edge  3  Rx 

Figure  D.l  Etotai  at  the  Rx  for  multiple  subpath  edges 


many  possible  paths  one  could  consider.  Hence,  a  new  method  that  derives  lit  diffraction 
path- loss  over  multiple  subpath  obstacles  using  path  clearance (/ij.)  and  the  first  Fresnel 
zone  radius(F'i)  is  simple  to  use. 
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D.2  A  Sample  Approximation  for  Lm 

The  field  intensity  beyond  a  subpath  obstacle  depends  on  the  form  of  the  obstacle. 
The  loss  due  to  a  knife-edge  obstacle  is  6  dB,  but  the  corresponding  value  for  a  smooth 
earth  is  about  20  dB  [Bullington,  1977],  Fig.D.2.a  shows  the  lit  diffraction  path-loss  for 
different  diffraction  parameter,  u,  where 

u  =  (D.3) 

and 


hf.  :  Path  clearance  (Fig4.2) 

Fi  :  Radius  of  First  Fresnel  zone. 

Fig.D.2.b  shows  the  path-loss  as  a  function  of  Note  that  the  lit  diffraction  path-loss 
linearly  decreases  as  ^  increases  until  ^  equals  0.6.  To  avoid  path-loss  a  clearance  of 
about  (0.6)Fi  or  more  is  required.  Thus  if  a  clearance  of  about  (0.6)jPi  or  more  exists, 
then  the  diffraction  effect  of  the  subpath  obstacle  can  be  ignored. 

Based  on  these  ideas,  I  suggest  the  following  equation  for  estimating  lit  diffraction 


Lit  Diffraction  Path-Loss  Lit  Diffraction  Path-Loss 


(a)  (b) 

Figure  D.2  Lit  diffraction  path-loss  for  different  diffraction  parameter 
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path-loss  due  to  obstacles  with  path  clearance(^c)  and  first  Fresnel  zone  radius(F'i): 


LiuidB) 

=  6-(^)10 

for  0  <  he  <  i0.6)Fi 

(D.4) 

LutidB) 

=  0 

for  he  >  (0.6)Fi 

(D.5) 

In  multiple  subpath  obstacle  situations,  only  the  obstacle  that  is  closest  to  the  LOS  (line  of 

sight)  projection  between  the  Tx  and  the  Rx  will  be  considered.  All  of  the  other  obstacles 

will  be  ignored.  Hence, 

I X  ■  dt  ■  dr 

(D.6) 

Fi 

y  dt  +  dr 

LutidB) 

=  6-(|r)10 

for  0  <  he  <  (0.6)F’i 

(D.7) 

where 

•  he'.  Height  difference  between  LOS(line  of  sight)  and  terrain  height  for  the  obstacle 
that  is  closest  to  the  LOS. 

•  dt'.  Distance  between  Tx  and  the  obstacle  that  closest  to  the  LOS. 

•  dr'.  Distance  between  obstacle  and  Rx. 

Fig.D.3  shows  the  construction  for  he,  dt,  and  dr  for  Eq.D.6  when  there  are  multiple 
subpath  obstacles. 


Terrain  Condition  between  Tx  and  Rx  One  Knife-Edge  between  Tx  and  Rx 


Figure  D.3  Find  lit  diffraction  knife-edge  from  real  terrain 
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D.3  Difference 


(a) 


(b) 


Figure  D.4  Etotai  at  the  Rx  by  the  different  consideration,  (a).  Consider  multiple  subpath 
obstacles  (b).  Ignore  the  obstacles  {Free  space) 


In  Fig.D.4,  if  there  are  muliple  subpath  edges  between  the  Tx  and  the  Rx  (Fig.D.4.a), 


then  Etotai  at  the  Rx  is 


{Etotai)  A  —  E-, 


incidence 


"b  {Edif fraction)  1  "b  {Edif fTaction)2  "b  {Edif fraction) Z 


(D.8) 


where 


{Etotai) A  •  Total  field  at  the  Rx  in  actual  communication 
{Ediffraction)n  ^  Diffracted  field  by  the  n-th  subpath  obstacle 

However,  if  the  effects  of  muliple  subpath  edges  are  ignored  (Fig.D.4.b),  then  Etotai 
at  the  Rx  is 

{Etotai)  F  Eincidence  (^•^) 

where 


{Etotai) F  :  Total  field  at  the  Rx  through  free  space 
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Finally,  using  Eq.D.8  and  Eq.D.9  ,the  Lm  is  estimated. 


{Etotal)F 


2 


(D.IO) 

(D.ll) 


The  Eq.D.ll  represents  the  basic  equation  which  estimates  over  multiple  subpath 
obstacles.  Fig.D.5  shows  difference  in  decibel  between  with  basic  theory  and  with 
suggested  equation.  The  difference  between  two  method  is  very  small  over  0  <  hc/Fi  < 
(0.6).  Hence,  through  research  suggested  equation  (Eq.D.7)  will  be  used  for  estimating 
L/ii  over  multiple  subpath  obstacles. 
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Appendix  E.  MATLAB  Codes 


The  MATLAB  codes  used  in  this  thesis  can  be  found  on  the  world  wide  web  at 

www.afit.af.mil/Schools/EN/ENG/LABS/C3EMRT/jyc.htm. 


1.  RTDMOD 

•  Estimates  propagation  path-loss  from  real  terrain  datas 

•  PL-total.m  :  For  n  by  1  input  data 

•  jangl.m  :  For  n  by  m  input  data 

2.  UTCMOD 

•  Estimates  propagation  path-loss  by  using  the  terrain  statistics 

•  PL-sto.m  :  For  n  by  1  input  data 

•  jang2.m  :  For  n  by  m  input  data 

3.  UTCMOD  for  SS 

•  Estimates  propagation  path-loss  of  spread  spectrum  communication 

•  PL-SS.m  :  For  n  by  1  input  data 

•  jangss.m  :  For  n  by  m  input  data 
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